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I.     INTRODUCTION 

A.     BACKGROUND 

Microwave  circuit  elements,  such  as  semiconductor  devices,  cannot  be  connected 
directly  to  the  terminals  of  a  microwave  network,  analyzer  for  measurement  of  scattering 
or  impedance  coefficients  but  must  be  coupled  through  a  chain  of  hardware  elements 
comprising  connector-adapters,  microstrip  sections,  and  bondwire  connections.  The 
basic  problem  under  consideration  is  the  removal  of  the  effects  of  intervening  hardware 
from  measurements  made  at  the  external  terminals  of  an  embedded  device  or  circuit 
component  to  obtain  the  intrinsic  device  characteristics.  Accurate  measurement  data 
on  a  device  is  essential  for  successful  circuit  design. 

Earlier  work  on  device  de-embedment  has  been  coupled  with  undesirable  restrictions 
such  as  the  assumption  of  physical  symmetry  of  hardware,  inadequate  provision  for  re- 
alistic measurement  standards,  or  neglect  of  statistical  treatment  of  data  for  improve- 
ment of  accuracy.  Consequently,  a  need  exists  for  a  practical  treatment  of  this 
measurement  problem  which  can  be  implemented  in  typical  laboratory  circumstances, 
with  acceptable  accuracy  of  results. 

In  order  to  measure  the  scattering  parameters  of  a  device,  the  device  is  often 
mounted  in  a  fixture.  This  fixture  has  a  set  of  scattering  parameters  associated  with  it 
so  that  the  parameters  measured  at  the  terminals  of  the  fixture  are  the  scattering  pa- 
rameters of  the  device  transformed  by  the  embedding  network  of  the  fixture. 

Thus,  to  obtain  the  scattering  parameters  of  the  device  itself,  the  device  must  be 
"de-embedded"  from  the  fixture.  Techniques  for  using  a  network  analyzer  (N.A.)  to  de- 
embed  the  device  from  the  test  fixture  have  been  discussed  earlier  (Refs.  1,2,10,11]. 

As  long  as  the  device  interface  connection  is  coaxial  line  or  waveguide,  de- 
embedding  is  relatively  straightforward  because  standards  are  available.  But  if  the  device 
is  constructed  without  a  coaxial  or  waveguide  interface,  as  in  the  case  of  microstrip  or 
coplanar  line  systems,  de-embedding  becomes  more  difficult. 

Extracting  the  device  parameters  from  the  overall  measured  parameters  is  strongly 
dependent  on  the  ability  to  characterize  the  fixturing.  Figure  1  is  the  block  diagram  and 
signal  flow  graph  of  a  typical  N.A.-  fixturing  setup.  The  problem  is  to  find  the  best 
characterization  networks  A  and  B  so  as  to  remove  their  contributions,  yielding  the  de- 
vice parameters,  S. 
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Figure  1.       Block  Diagram  and  Signal  Flow  Graph  of  a  Typical  Fixturing  Set-up 


B.    SCOPE 

The  scope  of  the  proposed  research  will  be  to  develop  a  working  procedure  for  car- 
rying out  the  removal  of  embedded  device  scattering  matrices  from  network  analyzer 
data.  The  procedure  developed  will  be  supported  by  experimental  measurements  on 
typical  microwave  circuits  fabricated  on  planar  substrates. 


II.     CHARACTERIZATION  OF  THE  FIXTURES 

A.     CHARACTERIZATION  CHOICES 

Several  papers  present  characterization  techniques  that  are  in  use  at  various  labo- 
ratories. There  are  four  main  choices  of  N.A.-fixturing  characterization  as  shown  Fig. 
2[Ref.  11. 


-    CHARACTERIZATION    CHOICES 

' 

i    Method 

1 

Method  2  , 

f 

.  Through 

(1).    Short  ckt 

S  -  parameters 

(2).  4  Open  ckt 

' 

fitted  to  a 
model 

1 

,   Method  3 

i    Method  4 

(1).  Through 

•  More  than  3 

(2).  Short 

known 

(3).  Delay 

refcoeff  of 

(Speciale) 

each  port 

Figure  2.      The  Four  Main  Choices  of  N.A.-Fixturing  Characterization. 


Method  (1)  combines  a  variety  of  measurements,  including  the  transmission  between 
fixtures  connected  back-to-back,  and  reflection  measurements  on  the  fixtures  terminated 
in  microstrip  stubs.  In  this  way,  data  is  assembled  to  enable  the  scattering  coefficients 
of  the  embedded  device  to  be  converted  to  the  device  scattering  matrix.  Method  (1) 
therefore  contains  features  of  both  methods  (3)  and  (4)  below,  and  represents  also  an 
early  example  of  the  use  of  microstrip  stubs  as  known  reference  standards. 

In  method  (2),  the  measured  S-parameters  are  fit  to  an  assumed  circuit  model  to 
calculate  the  device  parameters. 


Method  (3)  has  the  advantage  that  only  three  simple  test  pieces  in  the  device  me- 
dium are  required.    More  details  are  given  in  Chapter  V,  part  B. 

From  the  above  choices,  method  (4)  was  selected  for  investigation.  In  this  method, 
three  or  more  known  reflection  coefficients  are  measured  at  each  port.  The  reference 
load  reflection  coefficients  are  derived  from  independently  measured  or  calculated  non- 
zero values.  Reference  2  used  a  chip  capacitor,  wire  bonded  to  microstrip  on  alumina, 
with  capacitance  measured  at  low  frequency  (LF).  The  overall  reflection  coefficient  was 
then  measured  versus  frequency.  To  obtain  multiple  calibration  measurements,  a  piece 
was  scribed  off  the  chip,  and  the  measurement  cycle  repeated  twice  more  at  each  port. 

In  other  cases,  the  reference  load  was  a  high-Q  varactor  chip,  again  measured  as  a 
function  of  three  values  of  reverse  bias  at  LF.  The  chip  was  then  wire-bonded  to  the 
microstrip  conductor  and  the  overall  reflection  coefficients  of  the  fixture-chip  combina- 
tion versus  frequency  were  measured  at  each  value  of  reverse  bias. 

In  the  present  work,  the  reference  load  is  a  section  of  the  microstrip  line  itself.  Three 
different  lengths  of  microstrip  line  yield  the  three  different  reflection  coefficients  for  each 
port  by  separate  measurements. 

B.     THE  TEST  FIXTURE 

The  S-parameters  of  the  test  fixture  can  be  determined  by  measuring  a  set  of  known 
impedances  in  place  of  the  embedded  device.  This  process  is  known  as  "unterminating." 
Once  the  fixture's  S-parameters  have  been  determined,  the  device  is  measured  in  the 
calibrated  test  fixture,  and  the  device  parameters  are  "de-embedded"  from  the  data  to 
give  an  accurate  picture  of  device  performance. 

The  signal  flow  diagram  and  block  diagram  in  Fig.  3  define  the  input  and  output 
characteristics  of  the  test  fixture. 

Three  sets  of  S-parameters  must  be  determined  when  an  unknown  device  is  charac- 
terized as  a  two-port  network  mounted  in  a  test  fixture.  First,  the  S-parameters  of  the 
input  (A)  and  the  output  (B)  sides  of  the  fixture  must  be  evaluated  before  the  S- 
parameters  of  the  device  itself  can  be  deduced.  If  the  test  fixtures  are  known  to  be  sym- 
metrical, then  only  one  measurement  is  needed  to  determine  the  S-parameters  of  the  two 
sides  of  the  fixture. 

In  the  present  case,  however,  unsymmetrical  fixtures  are  chosen,  requiring  six  dif- 
ferent measurements  and  calculations  for  both  fixtures  independently.  To  solve  this 
problem,  the  input  half  of  the  fixture  is  first  terminated  at  the  device  plane  with  an  ar- 
bitrary load  of  reflection  coefficient,  fL,  followed  by  other  loads  with  different  reflection 
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Figure  3.      The  Flow  Diagram  and  Block  Diagram  of  the  Fixture. 


coefficients.    Then  the  output  half  of  the  fixture  must  also  be  considered  with  same 
method. 
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Figure  4.      One  Port  De- Embedment 

For  the  one  port  case,  as  shown  above,  the  reflection  coefficient  is  the  data  needed 
since  the  impedance  can  be  calculated  from  it: 


1  +rL 


(2.1) 


From  the  definition  of  the  scattering  matrix,  the  input  reflection  coefficient  is: 
a,5n  +a2Sn 


(2.2) 


where,  for  load  Tu 


a2  =  rLb2 


(2.3) 


and 


b2  =  S2xa\Jr  S22a2 


=  S2xax  +  S22{rLb2) 


(2.4) 


Thus, 


b2  =  S2l  '  (2.4a) 

Using  the  above,  we  get: 

Here,  Sn  can  be  found  by  measuring  a  matched  load  in  the  fixture  (TL  =  0).  If  we 
assume  that  this  measurement  gives  Tm  =  TML,  then  we  have  Sn  =  T ML.  Using  the  result 
in  Eq.  (2.5), 

(rln-rMl)(\-s22rl)^sns2lrL  (2.6) 

From  this, 

SnSifL  +  s22rL(rln  -  rML)  =  (r^  -  rML)  (2.7) 

Two  more  measurements  with  known  loads  will  provide  enough  data  to  determine 
S22  and  {SnS21).  Suppose  we  measure  with  an  open-circuited  load  (rL  =  +l).  and  a 
short-circuited  load  (rL  =  —1).  Then, 

sns2l  +  s22(roc-rMl)  =  roc-rML  {for    rL  =  +i)  (2.8) 

-sl2s2l-s22(rsc-rMI)  =  rsc-rML  {for    vL  =  -\)  (2.9) 

These  can  be  solved  for 

c        roc  +  Tsc  -2TML 

^22-  r  r  1>-1UJ 

1  oc~  l  sc 

c    c        ^   ^oc  ~  ^mD  (y  ml  "  r5C)  om 

■^12^21  — l  r  r  \z-il) 

1  oc~  l  sc 

With  these  two  quantities  known,  the  equation  for  de-embedding  a  general  unknown 
can  be  written  from  Eq.  (2.7) 


(r/„  -  r,VL) 

S12521  4-  522  (r,„  —  V ML) 


rL=  -  -    ,  %   /r:/L  -    ,  (2.12) 


But  in  the  actual  case,  it  is  very  difficult  to  construct  the  matched  load  and  ideal 
shorted  and  open  terminations  required  by  this  procedure.  Therefore,  we  rearrange  Eq. 
(2.12)  to  matrix  form  for  three  different  known  loads  and  their  reflection  coefficients. 
From  Eq.  (2.5), 
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1l 
r, 
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an  jii 


i  -  s22rL 


T       ^11        \w22  +  ^12^21 
1  L 

i  -  s22rL 


(2.13) 


Rearranging  this  gives: 


r, 


r       ~    r      +  */n^22  —  (^11^22  —  ^12^2l) 

1  L  l  L 


(2.14) 


—  +  r/ns22  -  a 


where, 


A  =  SUS22  —  SnS2\ 


(2.15) 


This  equation  then  involves  the  three  unknowns  5U,  S^,  and  A  .  These  three  unknowns 
can  be  determined  by  making  three  independent  measurements  of  r,„  with  three  different 
known  loads.    In  matrix  form,  the  problem  can  be  expressed  as: 


L\ 


r 


12 


U 


ln\ 


r,. 


in! 


r, 


M 


-1 
-1 
-1 


S22 
A 

= 

rL1 
rL2 
r^3 

(2.16) 


To  find  Su,  S22  and  A,  solve: 


u 


22 


1 

n«i 

-1 

rL, 

1 
rL2 

M*2 

-1 

i 

r» 

M/i3 

-1 

-1 

Mrtl 
W*2 

rL2 

Tl3 

(2-17) 


r^j,  TL2  and  rLi  are  three  different  calculated  reflection  coefficients  at  the  terminals 
of  the  three  different  loads.  Also,  Tml,  rin2  and  r,n3  are  three  reflection  coefficients 
measured  from  the  N.A.  with  respect  to  the  three  different  loads.  Now  the  problem  re- 
quires accurate  values  of  rLl,  ru,  Tu  and  TMl,  rM,  T^.  We  consider  in  the  next  chapter 
the  computation  of  ru,  rL2  and  TLi  for  reference  microstrip  line  loads. 
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III.     MICROSTRIP  LINE  FOR  THE  FIXTURE 

A.     MICROSTRIP  LINE 

The  definition  of  an  effective  relative  dielectric  constant,  ere,  allows  us  to  treat 
microstrip  line  like  a  TEM  system.  Relations  between  line  and  freespace  parameters 
then  hold  (which  are  typical  of  TEM  systems).  With  e„  replacing  the  conventional  rel- 
ative dielectric  constant, 


(3.1a) 
(3.16) 


S-ms 

1 

— 

, 

'■0 

V£r. 

V/715 

I 

c 

V  Cre 

Zo 

1 

^Oair 

v^7 

1 

vms   ^ms 

v    C 

(3.1c) 


(3.1<0 


where  Xms  and  vmt  are  the  wavelength  and  velocity  of  the  waves  on  the  microstrip  line  and 
/„  and  c  are  the  free  space  values  at  the  frequency  of  operation. 

In  Eqs.  (3.1),  Lms  and  Cms  are  the  inductance  and  capacitance  per  unit  length  of  the 
microstrip  line.  Z0  is  its  characteristic  impedance,  and  Z0air  is  the  value  Z0  would  have  if 
the  dielectric  layer  were  replaced  by  air.  To  the  lowest  order  of  approximation,  the  value 
of  er,  would  be  the  average  of  the  relative  dielectric  constant  of  the  substrate  and  that 
of  air,  namely: 

1  +e. 

lre  =  —^TL  (3-2) 


More  exact  values  have  been  obtained,  largely  through  solutions  of  the  quasistatic 
model,  corrected  to  match  measured  data. 
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Figure  5.      Microstrip  Line  Dimensions. 


w 


In  case  of  lines  with  line  width  ratio  u  =  —  =  0. 1  to  10, 

h 


£.  +  1        er  —  1 


LV 


1  + 


12.35 


+  0.06(1  -iOV0-44"*' 


(3.3) 


If  microstrip  were  a  true  TEM  system,  the  velocity  of  waves  on  the  line,  the 
impedance  Z0  and  other  line  parameters  would  be  independent  of  frequency.  The 
quasistatic  or  zero  frequency  values  of  the  parameters  are  good  approximations  at  fre- 
quencies up  to  the  order  of  a  few  gigaherz.  When  the  wavelength  becomes  comparable 
with  the  cross  sectional  dimensions  of  the  line,  v„,  and  the  line  parameters  are  no  longer 
frequency  independent.   Therefore,  dispersion  of  effective  relative  dielectric  constant  is 


treif)  =  *r 


h  ~  £re(°) 


fi,e(0)       2W    2 


(3.3a) 


where  h  is  the  substrate  height,  f  is  frequency  and  n0  is  the  permeability  of  free  space. 

B.     MICROSTRIP  STUB-RESONANCE  EVALUATION  OF  DIELECTRIC 
CONSTANT 

It  was  necessary  to  determine  the  substrate  dielectric  constant  e,  of  the  microstrip 
line  substrates  used  for  measurements  in  the  present  work. 
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The  steps  are: 

1.   Measure  the  resonant  frequency  of  a  quarter  wave  length  stub.  (Suhstrate  height 
is  h=  1.5  mm,  line  width  ratio  u  =  3.0  (where  u  =  w  /  h)), 


4  Different  microstrip  stubs,  with: 

—  Length  =  One  quarter  wave  length 

—  Substrate  height  =  1.5  mm 

—  w  /  h   (Line  width  ratio)  =  3.0 


connector 


|       L=Xn«/4 


T 


connector 


Figure  6.       Measurement  Setup  for  Dielectric  Constant  Measurement 

2.  Measure  the  value  of  52]  transmission  coefficient,  with  respect  to  frequency  using 
the  N.A., 

3.  Find  the  frequency  (f0)  at  which  maximum  dB  transmission  loss  due  to  resonance 
occurs,  as  shown  in  Fig.  7, 

S21 
$21"         Magnitude  plot 


Frequency 


Figure  7.       Determination  of  the  Max.  dB  Loss  Point. 
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Calculate  e,t  from  Eq.  (3.1a).  It  is  normally  required  to  account  for  the  fringing 
capacitance  at  an  open  microstrip  termination.  Computation  by  standard  methods 
[Ref.  3:  p.  138]  showed  that  the  open-end  length  correction  was  less  than  1.0%  for 
stub  lengths  greater  than  50  mm.  This  had  negligible  effect  on  the  calculated  £M 
and  was  neglected. 

;-0  =  }>msyJere  (3-4) 


\lLre 


therefore, 


;     = — 

where 


>-ms  =  4/- 


(3.5) 


(3.6) 


From  the  equations  above, 
Cre     [  4//0  ) 


(3.7) 


Through  the  several  measurements  as  shown  in  Table  1,  the  average  e„  ob- 
tained by  calculation  has  the  value  of  2.833., 


Table  1.     THE    MEASUREMENT    OF    EFFECTIVE    RELATIVE 
DIELECTRIC  CONSTANT 


Length  (mm) 

MGHz) 

£,. 

98.0 

0.452 

2.867 

80.0 

0.564 

2.763 

60.0 

0.746 

2.808 

40.0 

1.102 

2.895 

5.     From  stub  resonance  data,  the  corresponding  e„  is  calculated.    Eq.  (3.3)  was  in- 
verted by  trial  and  error  to  find  it  =  3.516  from  the  experimental  £,,. 

C.     CHARACTERISTIC  IMPEDANCE  OF  MICROSTRIP  LINE 

Expressions  for  the  characteristic  impedance  of  microstrip  lines  have  been  developed 
through  numerous  successive  approximations,  based  principally  on  TEM  models  having 
the  specific  microstrip  line  cross  section,  with  correction  terms  added  to  obtain  agree- 
ment with  experimental  observation. 
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v         60    , 

Z0  =  ^==  In 


(0.99u+ 1.795)(0.994w  +  5.S0)  1 


u{u+  1.268) 
where  the  synthesis  formula  for  linewidth  ratio  u  =  w/h  is: 


(3.8) 


"  =  f  =  TVUo.647  +  ^4U6.)+  1.321(1  +f)  (3.9) 


with 

._    v/2(fir+l)  1 

(3.10) 


A  = 


v'2(£.+  1) 


From  Eqs.  3.3,  3.6  and  3.8,  one  can  easily  get  the  characteristic  impedance  for  each 
microstrip  section.  For  example,  for  the  given  characteristic  impedance  of  50  Q  ,  the 
microstrip  line  width  ratio  is  u  =  2.283  and  the  width  is  3.43  mm.  In  the  construction 
of  microstrip  stub  loads,  it  is  necessary  to  choose  high  characteristic  impedance  for  ac- 
curate measurement,  since  high  impedance  leads  to  narrow  linewidth  and  a  well-defined 
point  of  attachment  of  the  stub.  Construction  problems  arise  in  making  microstrip  line 
with  very  small  width.  For  example,  if  we  want  to  make  a  microstrip  line  stub  of  90  Q 
characteristic  impedance,  the  required  width  is  about  1.2  mm  for  given  h  =  1.5  mm.  It 
is  difficult  to  handle  the  adhesive  foil  material  used  for  the  microstrip  itself.  For  this 
reason,  a  compromise  value  was  chosen,  using  2.0  mm  width  of  microstrip  stub  for  given 
substrate  height  h  =  1.5  mm.  Then, 

h  =  1.5mm 

w  =  2.0mm 

also  u  =  -r-=  1-33. 
h 

Therefore, 

Z0^  =  68.2Q, 

which  was  chosen  as  a  working  value  for  stub  impedance. 
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MICROSTRIP  UNE 

CHARACTERISTIC  IMPEDANCE  :  50  Ohm 
WIDTH  :  3.43  mm 

OPEN  STUB 


CHARACTERISTIC  IMPEDANCE  :  68.2  Ohm 
WIDTH :  2.00  mm 


connector 


68.2  ohm 


.2.0  mm 


Figure  8.       Experimental  Setup  for  Characterization  of  the  Fixture 

D.     STUB  LOAD  IMPEDANCE  AND  REFLECTION  COEFFICIENT 

Consider  the  transmission-line  impedance-transformation  law: 


£-ln  ~  M)j 


ZL  cos  (31  +jZ0S  sin  (31 
Z0s  cos  {31  +  ZL  sin  (31 


(3.11) 


From  Eq.  (3.11),  when  ZL  goes  to  oo    (open  ended),  the  input  impedance  of  an  open  - 
ended  stub  becomes: 
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Z/„   —  Z, 


05 


cos  PI 
_  j  sin  pi  _ 


(3.1. 


=  -JZ0s  cot(pf) 


using 


(3.13) 


/ 


(3.14) 


Therefore, 
ZL  =  -JZ0s  cot 


InfljT, 


=  -jz< 


Qs 


r  ^gvfrs  i 

C0SL~^7J 


(3.15) 


Once  the  Z0j  is  chosen,  we  can  determine  the  load  impedance  of  three  different  loads. 
For  example,  given  Z0,  =  68.2  ohms  and  microstrip  line  lengths  of  4  cm,  3  cm,  2  cm,  the 
load  impedances  are  tabulated  as  shown  Table  2. 

Table  2.      THE  LOAD  IMPEDANCE  OF  REFERENCE  STUB  LOADS 


Frequency(GHz) 

L  =  4  cm  (Z„) 

L  =  3  cm  (Zu) 

L  =  2  cm  (Zi3) 

1.00 

-j  11.06 

-j  38.44 

-j  80.15 

2.00 

j  204.8 

j  41.28 

-j  11.06 

3.00 

-j  35.69 

-j  1230 

j  41.28 

4.00 

j  91.05 

-j  35.69 

j  204.81 

5.00 

-j  70.73 

j  44.23 

-j  168.99 
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The  three  different  reflection  coefficients  of  the  given  arbitarv  loads  are: 


r„  = 


z, 


Z-Li  +  ^0 


(3.16) 


where  i  =   1,  2,  3.    Therefore,  from  the  given  stub  impedances  and  Z0  =  50Q,  the  re- 
flection coefficients  are  tabulated  in  Table  3. 


Table  3.      THE  STUB  REFLECTION  COEFFICIENTS 

Freq.(GHz) 

L  =  4  cm  (ru) 

L  =  3  cm  (rL2) 

L  =  2  cm  (TJ 

1.00 

-0.91  -J0.42 

-0.26  -  J0.97 

0.44  -  J0.90 

2.00 

0.89  +  J0.46 

-0.19-Hj0.98 

-0.91  -  jO.42 

3.00 

-0.32  -  J0.95 

1.00- jO.05 

-0.19-t-j0.98 

4.00 

0.54  +  J0.84 

-0.32  -  J0.95 

0.89  +  J0.46 

5.00 

0.33  -  jO.94 

-0.12-Hj0.99 

0.84  -  J0.54 
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IV.     APPLICATION  OF  THE  DE-EMBEDDING  PROCEDURE 

A.  MEASUREMENT  OF  THE  EMBEDDING  FIXTURE 

The  reflection  coefficients  required  are  measured  by  the  HP  8409B  vector  network- 
analyzer.  One  can  find  the  reflection  coefficient,  rm,  of  the  stub  terminated  fixtures  by 
measurement  of  the  scattering  matrix  element  Sn  using  the  network-analyzer.  The  out- 
puts of  the  N.A.  represent  the  magnitude  (dB  loss)  and  phase  angle  (degrees)  of  Sw 
Therefore,  the  outputs  of  the  N.A.  must  be  changed  to  complex  notation  using  the 
simple  modification: 

~SUJB  ~Sl\dB 

r//7  =  (10)    20     cos(argSn)+J(10)    20     sin(argSn).  (4.1) 

The  outputs  of  the  N.A.  and  TM  are  obtained  for  each  load  on  one  fixture  (A)  as  shown 
in  Figs.  9-a  to  9-c.  When  we  consider  the  other  fixture  (B),  all  steps  are  exactly  same 
as  for  fixture  A  except  the  fixture  shape  itself,  shown  in  Figs.  10-a  to  10-c.  Tables  for 
fixture  A  and  fixture  B  data  are  available  in  Appendix  A. 

B.  S-PARAMETERS  OF  THE  FIXTURE 

To  consider  the  de-embedding  calculation,  we  return  to  the  original  matrix  equation. 
The  success  of  this  characterization  technique  depends  heavily  on  the  measurement  of 
the  load  impedances  used  in  the  unterminating  process  for  the  fixtures.  Once  an  input 
reflection  coefficient  has  been  measured  with  three  different  impedance  standards,  the 
S-parameters  of  the  test  fixture  can  be  computed  using  the  matrix  introduced  earlier. 
A  word  of  caution  must  be  offered  about  the  convergence  of  this  apparently  simple 
problem. 

The  matrix  that  must  be  inverted  consists  of  elements  very  near  unity,  a  situation 
often  used  to  illustrate  some  of  the  numerical  inversion  methods.  In  this  case,  the 
Gauss-Jordan  method  with  full  pivoting  is  recommended.  Pivoting  improves  the  stabil- 
ity of  the  numerical  solution  by  interchanging  rows  and  columns  in  the  Gaussian  elim- 
ination. In  the  notation  of  Eq.  (2.17),  the  scattering  parameter  S,y  is  changed  to  AtJ  and 
BtJ  for  fixture  A  and  B  (respectively): 
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Figure  9.       Measurement  of  the  Load  for  Fixture  A 
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Figure  10.       Measurement  of  the  Load  for  Fixture  B 
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(4.2) 


where 


A  =  AUA22  —  AnAlx. 


(4.3) 


The  computed  values  for  An  and  An,  combined  with  a  knowledge  of  A  ,  can  be  used  to 
determine  the  An  A2l  product.  Note  that  the  fixture  is  reciprocal,  so  Al2  =  An  and  for  the 
other  fixture,  Bl2  =  B21.  Thus,  the  magnitudes  of  these  terms  can  be  calculated  by  taking 
the  square  root  of  the  product  term  in  the  expression  for  A  . 

The  question  of  phase  is  left  open,  but  is  not  essential.  Since  the  Al2,  A2U  B12  and 
B2l  terms  always  appear  in  product  form,  it  does  not  matter  which  of  the  two  possible 
phase  angles  one  chooses,  as  long  as  the  phase  angles  for  A2X  and  B2l  are  approximately 
the  same.  All  of  these  calculations  are  shown  in  Appendix  E  using  the  Fortran  program 
shown.  This  program  combines  the  measurement  and  calculated  values,  that  is,  the 
outputs  of  the  N.A.  are  used  in  the  calculation  as  a  form  of  datafile.  The  resulting  S- 
parameters  of  these  fixtures  are  shown  in  Appendix  B. 

C.     DE-EMBEDDING  PROCEDURE 
1.     Background 

When  a  two-port  network  is  connected  between  fixture  A  and  B,  the  signal  flow 
graph  for  the  system  may  be  drawn  as  shown  in  Fig.  11.  Parameters  AB  and  BA  are 
mainly  contributed  by  the  finite  isolation  of  the  internal  switches  of  the  network  ana- 
lyzer. The  terms  Su,  S2U  Sl2,  and  S22  are  the  scattering  parameters  of  the  device  under 
test.  The  fixture  parameters,  An,  A2l,  Al2  and  Au,  can  be  determined  by  reflection 
measurements  at  port  1,  as  discussed  earlier.  Parameters  B22,  Bn  and  BnB2{  can  be  eval- 
uated by  similar  measurements  at  port  2.  Two  port  transmission  measurements  with  the 
device  under  test  replaced  by  a  through-line  are  used  for  evaluating  parameters  AB  and 
BA,  and  products  A2lB2l  and  AnBx2.  When  a  through-line  of  length  /,  having  zero  re- 
flections at  both  ends  and  a  transmission  T  =  e~JI",  is  connected  between  the  two  ports, 
the  signal-flow  graph  of  Fig.  11  can  be  simplified  to  the  one  shown  in  Fig.  12. 
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Figure  11.       Signal  Flow  Graph  for  Measurements  on  a  Two  -  Port  Network 


The  measured  transmission  Tm  can  be  expressed  as 


Tm  =  BA  + 


\-A22BuT2 


(4.4) 


When  the  error  parameters  A22  and  Bn  have  been  determined  by  the  reflection  measure- 
ments, two  transmission  measurements  must  be  performed  for  the  determination  of  BA 
and  A2lB2l.  The  value  of  BA  can  be  found  by  making  a  transmission  measurement  when 
no  connection  exists  between  the  ports  (that  is  T  =  0).  The  parameter  A2lB2l  can  be 
determined  from  the  measurement  of  a  through  -  line  using 


A2lB2l  = 


{Tm-BA){\-A22Bue-npl) 


(4.5) 


where  T=e~'fil,  j/?  being  the  propagation  constant  of  the  through-line.    The  error  pa- 
rameters AB  and  Al2Bl2  can  be  determined  by  making  a  similar  measurement. 

Now  that  the  scattering  characteristics  of  the  embedding  networks  have  been 
measured,  the  scattering  parameters  Afllf  Mu,  Ml2  and  M2l  can  be  corrected  for  the  em- 
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Figure  12.       Signal  Flow  Graph  with  a  Through  -  Line 

bedding  networks  to  obtain  Sn,  Sn,  Sl2  and  S2l.  Elements  Mtj  represent  the  scattering 
parameters  of  both  fixtures  and  device,  namely: 


\fn  =  An  + 


^1 1  ^  12^2l( ^  ~  ^22^1l)  +  (S2\S\2S\\A\2A2\) 

D 


(4.6) 


Jiyf  D   <     ,      ^21^21^21 

M2l  =  BA+ ^ 


(4.7) 


M22  =  B22  + 


Mn=AB+ 


S22B2[B]2(l  -  £11^22)  +  (^21^12^22^12^21) 
D 


SnBnAn 
D 


(4.8) 
(4.9) 


where 


D  =  \  -  Sn^22  —  ^22^11  -  ^12^21^22^11  +  ^22^11^22^11 


(4.10) 
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These  equations  must  be  inverted  to  extract  the  scattering  parameters  S,r  An  explicit 
solution  for  the  scattering  parameters  of  the  device  under  test  has  been  presented  in  Ref. 
[4]  :  p.  313  .   The  solution  is  as  follows: 

G(\+BUH)-BUEF 
Sn  = "y (4.11) 

Si2  =  —  (4.12) 

S2i=—  (4.13) 

_         H([  +  A22G)  -  A22EF 

522= Jj (4.14) 

where 

N=  (1  +  A22G)(l  +  BUH)  -  AnBuEF  (4.15) 

(Mn-AB) 

(Af2,  -  5/1) 

G=_7^r  (418) 

(A/22  -  fl22) 

H-uhr\  (4-19) 

In  the  present  measurement,  inclusion  of  the  leakage  paths  BA  and  AB  is  not  necessary 
because  of  the  one-port  calibration  procedure  used.  Therefore,  Eqs.  (4.16,  4.17)  are: 


(4.16a) 

(4.17a) 
Some  variations  of  the  above  calibration  procedure  are  available  in  Ref.  [4]  :  p.  315  . 
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E  = 

Mn 

^12^12 

n 

M2l 

t  — 

^21^21 

The  de-embedding  procedure  used  in  this  work  is  based  on  Eqs.  (4.11  to  4.15) 
and  Eqs.  (4.16a,  4.17a,  4.18  and  4.19)  which  use  the  9  key  parameters: 

S\\,  S{2,  S2\,  S22,  E,  F,  0,  II,  N 

and  the  foregoing  derivations.  The  complete  block  diagram  is  shown  in  Tig.  14.  If  we 
know  A,j  and  Btl,  the  scattering  parameters  of  fixtures  A  and  B,  and  the  parameters  of 
both  fixtures  plus  device,  MtJ,  then  the  scattering  parameters  of  the  device  itself  are  ob- 
tained by  using  Eqs.  4.11  to  4.14.  In  order  to  test  the  de-embedding  procedure,  consider 
the  set  up  of  Fig.  13. 
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FIXTURE  +  DEVICE 


FIXTURE  A 


CONNECTOR 


N\\\\\\\\\\\\\\\\\\\\N 


\ 


MICROSTRIP  LINE 


Figure  13.      Device  and  Fixture  Assembly  for  the  51  Ohm  Resistor 

It  is  necessary  to  choose  a  device  that  has  known  scattering  parameters;  then  it  will  be 
possible  to  compare  the  output  of  the  de-embedding  experiment  with  the  theoretical 
parameters  of  device  itself.  The  initial  test  object  is  a  series-connected  51  Ohm  resistor 
which  has  a  small  inductance  so  that  the  impedance  of  the  51  Ohm  device  is: 


Z=  R+jcoL 
w  51  +jlnfL 

where  we  may  assume  the  parasitic  inductance  associated  with  the  resistor  to  be 


(4.20) 
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Figure  14.       The  Block  Diagram  of  the  De- Embedding  Procedure 
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L  =  2.5  nH. 

The  scattering  matrix  of  an  impedance  Z  connected  in  series  with  the  trans- 
mission line  is  : 


5    = 


^11         ^12 
S2l        S22 


2  +  z 

2 


2  +  z         2+z 


(4.21) 


where 


Zn 


and  Z0  is  the  characteristic  impedance  (50  Ohm)  of  the  connected  line.   In  this  case,  the 
TOUCHSTONE  CAD  program  is  useful  in  calculating  the  scattering  parameters. 

These  theoretical  outputs  are  shown  in  Appendix  C.    Also  the  magnitude  plots 
of  the  calculated  scattering  parameters  of  this  device  are  shown  in  Fig.  15. 
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Figure  15.       Theoretical  Scattering  Parameters  of  51  Ohm  Resistor 
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2.     Device  Embedment 

The  S-parameters  of  both  fixtures  plus  device  are  measured  by  the  network  an- 
alyzer at  planes  AA'  and  BB'  as  shown  in  Fig.  16.  Typical  measured  magnitude  data  is 
shown  in  Fig.  17  with  data  tables  in  Appendix  C. 
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Figure  16.      Signal  Flow  graph  for  the  Measurement  of  the  (Fixture  +  Device)  S- 
parameters 
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Figure  17.      The  Magnitude  Plot  of  Both  Microstrip  Fixture  and  51  Ohm  Resistor 
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3.     De-embedding  the  Device 

Now  we  are  ready  for  calculation  of  the  S-parameters  of  the  device.  The  S- 
parameters  of  fixture  A,  fixture  B  and  both  fixtures  plus  device  constitute  the  necessary 
data.  Once  the  microstrip  circuit  board  is  prepared,  then  the  network  analyzer  provides 
the  three  different  input  reflection  coefficients  r,„„  rM  and  Tm3  with  respect  to  three  dif- 
ferent calibration  loads.  Also  computations  are  carried  out  to  find  the  load  reflection 
coefficients  of  these  three  different  loads,  Tu,  rL2  and  Tu.  After  following  the  same 
procedure  for  both  fixture  A  and  fixture  B,  we  can  get  the  respective  S-parameters. 
From  the  matrix  equation  (Eq.  2.17),  the  three  unknowns  Sn,  S22  and  A  are  solved.  For 
de-embedding  the  device  itself,  we  need  the  additional  S-parameters  for  the  device  plus 
Fixture.  These  three  sets  of  S-parameters  are  utilized  to  solve  the  de-embedding  problem. 

Equations  4.11  to  4.14  are  now  ready  for  use.  Through  these  equations,  the 
de-embedded  outputs  are  obtained,  as  shown  in  Fig.  18  (the  table  of  de-embedded  S- 
parameters  is  shown  in  Appendix  C).  In  addition  to  the  plot  of  the  de-embedded  mag- 
nitudes, the  comparison  of  the  scattering  parameters  Sn,  Sl2  and  Su  of  the  de-embedded 
output  with  the  parameters  of  the  embedded  resistor  and  the  theoretical  values  are  also 
available  in  Fig.  19  to  Fig.  21  (data  tables  in  Appendix  C). 
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Figure  18.       Magnitude  Plot  of  the  De-embedded  Parameters  of  51  Ohm  Resister 
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Figure  19.      Comparison  of  the  Sll  of  the  51  Ohm  Resistor  with  the  Theoretical 
Values 
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Figure  20.       Comparison  of  the  S12  of  the  51  Ohm  Resistor  with  the  Theoretical 
Values 
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Figure  21.      Comparison  of  the  S22  of  the  51  Ohm  Resistor  with  the  Theoretical 
Values 
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V.     FURTHER  DE  -  EMBEDDING  METHODS 

A.     MEASUREMENT  OF  DEVICE  PARAMETERS  USING  A  SYMMETRIC 
FIXTURE. 

1.     Background 

As  discussed  earlier,  if  the  device  is  constructed  without  a  coaxial  or  waveguide 
interface,  de-embedding  becomes  more  difficult.  In  some  cases,  it  is  possible  to  proceed 
using  an  auxiliary  standard  whose  properties  are  calculated  or  measured  at  low  fre- 
quencies, but,  in  other  cases,  no  suitable  auxiliary  standard  is  available.  In  either  case, 
if  the  fixture  is  symmetric,  a  de-embedding  procedure  that  does  not  require  the  use  o^ 
an  auxiliary  standard  can  be  developed  [Ref.  5]. 

The  concept  for  symmetric  de-embedding  is  illustrated  in  Fig.  22.  Fig.  22-a 
shows  a  device  inserted  in  a  symmetric  test  fixture.  We  assume  that  the  device  has 
complementary  interfaces  on  the  input  and  output  sides  so  that  when  the  device  is  re- 
moved from  the  fixture,  the  two  halves  of  the  fixture  can  be  connected,  as  shown  in  Fig. 
22-b.  The  scattering  parameters  of  each  half  of  the  test  fixture  can  be  extracted  from  the 
scattering  parameters  of  the  whole  test  fixture  connected,  as  illustrated  in  Fig.  22-b. 
Then  the  device  is  inserted  into  the  test  fixture,  as  illustrated  in  Fig.  22-a,  and  the  scat- 
tering parameters  are  again  measured. 

The  method  for  implementing  the  concept  outlined  above  is  as  follows.  If  we 
have  a  fixture  whose  T-parameters  in  a  50  ohm  system  are  ZTF~],  then  if  [_HTF~\  re- 
presents the  square  root  of  [_TF~\, 

lHTr\LHTFl  =  lTr\.  (5.1) 

Now  suppose  we  want  to  measure  the  T-parameters  of  a  device  in  a  50  ohm  system 
CrSJ.  If  we  separate  the  fixture  into  two  symmetric  parts  and  inserted  the  device  to  be 
measured,  then  for  the  combination  we  would  measure  a  set  of  T-parameters  [TFS2 

given  by: 

lTFSl  =  iHTr\lTSHHTFl-  (5-2) 

To  find  CTS],  we  would  simply  invert  above  equation  to  get 

CTS]  =  ZHTFl~lZTFSlZHTFrl>  (5-3) 
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Figure  22.       Device  Embedded  Between  Symmetric  Fixture  and  Test  Fixture 

This  procedure  de-embeds  the  device  parameters  \TS~\  from  the  measured  parameters 
C  TFS~\  and  Q  TF~\.  The  steps  required  to  find  [  TS}  in  a  symmetric  fixture  are  as  follows: 

1.  Measure  the  scattering  parameters  of  both  fixture  plus  device,  CSFS],  using  the 
network  analyzer, 


tSFS} = 


SFSU 

SFS2l 


SFSn 
SFS22 


(5.4) 


then  transform  to  the  T-parameters: 
-A/,  Si 


ZTFS2  = 


'n 


SFS 


21 


-SFS 


22 


SFS 


21 


SFS2, 

1 
SFS2] 


(5.5) 


where  A/,  =  SFSnSFS22  +  SFSl2SFS2l. 

2.  Measure  the  scattering  parameters  of  both  fixtures  ZSF2.  In  practice,  both  fixtures 
are  measured  with  the  network  analyzer  individually: 
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ZHSF1  =  ZHSFll  =  IHSF22  = 


HSFU 
HSF2l 


HSF 


12 


HSF 


22 


Then,  transforming  to  the  T-parameters: 
-M7  HSF, 


tHTFJ  = 


HSF 


21 


HSF 


21 


-HSFn 


1 


HSF 


21 


//Si7- 


21 


where  M2  =  HSFnHSFa  +  HSFl2HSFu. 
3.   Invert  ZHTF3: 


£HFF}~]  = 


HTF 


22 


-//7F- 


21 


M3 
-HTFn 


M, 


//rf, 


A/, 


AA 


where  A/3  =  HTFnHFFn  -  HTFnHTF{ 
4.  Find  device  T-parameters  r_7Sl: 

1 


crs:  = 


HTFU    HTFn 


HTF,,    HTF 


21 


22 


TFSU    TFSn 


TFS2l    TFS 


22 


HTF,    HTF 


ii 


12 


-1 


HTF2X    HTF22 


5.  Transform  fTS]  to  S-parameters  [S51: 


tssi  = 

where 


S 


21 


>22 


s„  = 


TS 


12 


'11 


75 


22 


S12  —  TSU  — 

s L_ 


TSl2TS2l 


TS 


22 


s22  —  — 


22 

TS_ 
TS 


(5.6) 


(5.7) 


(5.8) 


(5-9) 


(5.10) 

(5.11) 
(5-12) 

(5-13) 
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2.     De-Embedding  of  a  Filter  Unit 

As  a  test  demonstration  of  the  symmetrical  de-embedment  procedure,  an  8436A 
bandpass  filter  (8-12  GHz)  is  used  as  the  test  device.  All  components  of  the  network 
had  type  APC-7  test  port  connectors.  The  outputs  CSFS]  and  [HSF2  of  the  network 
analyzer  are  shown  in  Appendix  D.  To  get  the  de-embedded  outputs,  the  Fortran  pro- 
gram shown  in  Appendix  F  is  used.  The  magnitude  plot  of  the  filter  embedded,  the  filter 
itself,  and  the  de-embedded  filter  are  shown  in  Figs.  23  to  25.  Data  Tables  are  shown 
in  Appendix  D.  The  comparison  of  the  scattering  parameters  Su,  S12  and  S22  of  the  filter 
itself,  and  embedded  filter  are  shown  in  Figs.  26  to  28.  Tabulated  data  is  shown  in  Ap- 
pendix D. 
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Figure  23.       The  S-paranieter  Magnitude  Plot  of  Both  Fixture  and  Filter 
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Figure  24.       The  S-parameter  Magnitude  Plot  of  Filter  Without  Fixtures 


dB 

Loss 


50 


40   • 


30 


20 


10 


\ 

\ 
\ 

De-embedded  Outputs 

V 

\ 

\ 

1    1 

Magnitude  plots 

• 

1 

\ 
\ 

\ 
\ 

T 

v 

H1    /^A 

-^ 

■^J  u 

S12    and    S21               „ , 

6.00      6.64      7.28      7.92      8.56      9.20      9.84      10.5      11.1        11.8      12.4 

Freq.  (GHz) 


Figure  25.       The  S-parameter  Magnitude  Plot  of  De-embedded  Parameters  of  Filter 


43 


dB 
Loss 


50 


40 


30 


20  ■ 


10 


Comparison  of  Output  S11 
(1).  Fixtures  +  Filter 
(2).  Filter  without  Fixtures 
(3).  De-embedded  output 


t 

/' 

4 

\ 
\ 

i 

i 
i 

\ 

\ 

i  . 

\ 

4-'. 

>Cs 


/ 


£*-: 


V 


v/ 


/ 


6.00      6.64      7.28      7.92      8.56      9.20      9.84      10.5      11.1        11.8      12.4 

Freq.  (GHz) 


Figure  26.      Comparison  of  the  Sll  of  De-embedded  Filter  »ith  Filter  Without  Fix- 
tures 
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Tigure  27.      Comparison  of  the  S 12  of  De-embedded  Filter  uith  Filter  Without  Fix- 
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Figure  28.      Comparison  of  the  S22  of  De-embedded  Filter  with  Filter  Without  Fix- 
tures 
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B.     TSD  (THROUGH,  SHORT,  DELAY)  METHOD 

This  procedure,  developed  by  R.  Speciale  [Rcfs.  6,  7  ],  is  intended  for  system  calibration 
and  systematic  error  removal  in  automated  S-parameter  measurement.  This  procedure 
uses  only  three  non-critical  reference  standards,  called  the  "Through",  the  "Short"  and 
the  "Delay",  and  leads  to  a  set  of  equations  for  which  a  fast,  non-iterative,  explicit  sol- 
ution is  known  [Ref.  6]. 
1.     The  Basic  Concept 

This  procedure  uses  only  three  mechanically  simple  reference  standards,  known 
respectively  as  the  "Through",  the  "Short",  and  the  "Delay". 

The  "Through"  is  defined  as  the  direct  mutual  connection  of  fixtures  A  and  B, 
as  shown  in  Fig.  29,  when  the  uncalibrated  S-parameter  data  for  the  device  are  acquired. 
The  S-parameter  data  for  "Through"  are  therefore  acquired  with  the  measurement  lines 
to  the  device  directly  interconnected  and  the  device  removed  from  the  signal  path. 

The  "Short"  is  defined  as  a  reference  object  presenting  perfect  shorts  at  the  im- 
mediate interfaces  where  the  unknown  device  is  connected  for  measurement. 

Finally,  the  "Delay"  is  defined  as  any  arbitrary  length  of  50  ohm  line  inserted 
between  the  interfaces  where  the  unknown  device  is  connected  for  measurement.  Sum- 
marizing, for  the  "Through"  and  the  "Delay",  the  signal  path  is  continuous  between  ports 
1  and  2  of  the  S-parameter  unit,  through  the  couplers  used  for  connecting  the  object. 
A  direct  connection  represents  the  "Through"  and  a  finite  arbitrary  length  line  for  the 
"Delay",  in  place  of  the  unknown  object.  On  the  contrary  the  signal  path  is  blocked  for 
the  "Short"  by  perfect  shorts  appearing  at  the  interfaces  where  the  unknown  device  is 
normally  connected. 

Three  S-parameter  matrices  C-5r"J,  [5S]  and  ZSD1  are  acquired  at  each  desired 
frequency  during  system  calibration.  These  matrices  describe  the  "Through",  the  "Short" 
and  the  "Delay"  respectively.  Because  of  the  device  imperfections,  these  matrices  do  not 
represent  theoretically  perfect  "Through",  "Short"  and  "Delay"  elements,  which  are  as- 
sumed to  exist  in  the  actual  measurement. 

Now  consider  the  error  boxes  A  and  B.  Error  box  A  is  assumed  to  be  connected 
on  the  side  of  port  1,  while  error  box  B  is  assumed  to  be  on  the  side  of  port  2,  as  shown 
in  Fie.  29. 
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Figure  29. 


Block  Diagram  of  the  TSD  Calibration  Procedure 


From  Fig.  29,  the  error  boxes  A  and  B  are  always  assumed  to  be  inserted  at  the 
two  ends  of  the  measurement  50  ohm  lines  connecting  the  device  to  ports  1  and  2  of  the 
S-parameter  unit.  An  algorithm  has  been  derived  wh  .h  performs  an  explicit,  non- 
iterative  calculation  of  the  S-parameter  matrices  ZSA3  and  CS^D,  which  characterize  the 
virtual  error  boxes  A  and  B. 

2.     The  Equations  and  their  Explicit  solution 

As  discussed  earlier,  the  T-matrices  can  be  multiplied  simply  with  cascaded  ele- 
ments.  Therefore,  consider  Fig.  30,  where  the  T-matrix  for  the  connection  shown  is: 


ITADB1-ZTA1  LTD1  LTB1. 


(5.14) 


From  the  above,  we  get: 


LTDl  =  LTAT]  LTADB1  crfi] 


-i 


(5.15) 
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Device  Embedded  Setup 


where  [  TD~\  is  the  corrected  T-matrix  and  [  TADB~\  is  the  uncalibrated  T-matrix  for  the 
device. 

The  T-parameters  of  fixtures  plus  device,  [_TADB],  are  given  by  the  network 
analyzer.  But  the  problem  is  to  obtain  the  T-parameters  of  each  fixture,  \_TA~\  and 
r_rB].  To  do  this,  we  consider  the  error  boxes  A  and  B  as  shown  in  Fig.  31. 

From  a  formal  point  of  view,  both  the  "Through"  and  the  "Delay"  can  be  con- 
sidered as  different,  finite  lengths  of  50  ohm  lines,  cascaded  between  the  error  boxes  A 
and  B.  With  \_TDl~\  and  [_TD2~\,  the  T-matrix  derived  from  the  S-parameter  data  obtained 
in  these  two  calibrations,  and  [7^3  and  LTB3,  the  unknown  T-parameters  of  the  error 
boxes  A  and  B,  we  can  write: 


ZTDll-ZTAlZTLlllTBl, 


(5.16) 


crn2i  =  :r.][r,2][rfl] 


(5.17) 


For  the  50  ohm  lines  L,  and  L2  we  have: 


ZTLJ  = 


,-p\ 


0 

oP\ 


(5.18) 


where  P\  =  h+jQx 
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Figure  31.       Error  Boxes  of  A  and  B  with  Delay 
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where  p2  =  r2  +jd2. 

Equations  (5.16)  and  (5.17)  together  with  the  two  equations  for  the  "Short" 
characterize  boxes  A  and  B.  Also,  consider  the  situation  for  the  "Short",  as  shown  in 
Fig.  32,  we  have 


•Ji  i  ^  — 


ns~  a, 


Si  15  _ 


b. 


IB 


11B 


SUA  4-  det  SA 

i  +  s22A 

S22B  +  det  SB 

l  +  S 


(5.20) 


(5.21) 


ns 


where  Sns  and  S^  are  Sn  and  Su  of  the    Short",  and  SilA,  SliB  are  Sa  of  box  A  and  box 
B  (det  means  matrix  determinant),  where  i=  1.  2. 
From  Eqs.  (5.16)  and  (5.17),  we  get: 


■->!  1  J    — 


n\A 


Tn 


(5.22) 


HA 
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Figure  32.       Short  Situation  of  Two  Error  Boxes 


det  Sa  =  5-,-, 


2  A 


Si 


(5.22a) 


JlIK  — 


—  Ti\b       —  £l 


lllB  K\ 


det  SR  =  S 


-a; 


B  ~  °  1 1 B        £ 


Now  substitute  Eqs.  (5.22)  and  (5.22a)  into  (5.20)  to  obtain: 

#3 


■^22.4  — 


#, 


~~  ^115 


•^HS  — 


Si 


;5.23) 


(5.23a) 


(5.24) 


and  also, 

$\2A$2\A  =  S\\A$22A  ~  ^et  $A 


■Jll   4 


'MA 


Si 


n2A 


(5.25) 
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Substituting  Eqs.  (5.23)  and  (5.23a)  into  (5.21), 


^4 


-s22 


225 


iiifl  — 


Sl2S  + 


A', 


(5.26) 


^\2B^2\B  ~  S\\B$22B  ~  ^et  $1 


$22B  +     £' 


s 


11B 


(5.27) 


where 


*'-* 

Hn-Hu 

I  H — 

2  sinh  p 

«--± 

K22  -  Kx  | 
2  sinh  p 

//12 

2  sinh  p 

A'3  =  -- 

•M2 

J  sinh  p 

ff4  =  -- 

H2] 

2  sinh  p 

is                    -1 

4   '  2  sinh 


(5.28) 


(5.29) 


(5.30) 


(5.31) 


(5.32) 


(5.33) 


Here  the  quantities  H,,  and  Kv  are  the  elements  of  the  matrices  C#D  and  [/Q,  defined 
bv 


c//]  =  :rflj  ITD21 


-i 


(5.34) 


[Ai  =  cr02r1  crD1: 


(5.35) 


and  p   is  defined  by: 
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=  cosh  1[y(^n  +  ^22)]  = 


T(KU  +  K 


a] 


(5.36) 


Once  the  S-parameters  of  the  error  boxes  A  and  B  are  known  as  given  by  Eqs.  (5.22) 
through  (5.27),  use  the  general  conversion  formulas, 


:n  = 


T\\      Tn 


^21         ^22 


Si 


-detS      S, 


S 


22 


I 


(5.37) 


where  det  5  =  SnS22  —  Sl2S2i. 

As  discussed  before,  the  following  calculation  to  de-embed  the  device  matrix, 

LTD1  =  LTAT[     LTADB1     ITBT\ 

can  be  carried  out  using  Eqs.  (5.22)  through  (5.27). 
3.     Discussion 

This  calibration  procedure  uses  only  three  reference  standards,  which  is  the 
minimum  necessary  for  characterizing  the  error  boxes  A  and  B.  By  using  two  delays  of 
different  length,  for  a  total  of  4  standards,  a  check  of  intrinsic  consistency  may  be  ob- 
tained. The  standards  used  are  simple  enough  to  be  inserted  at  the  location  of  the  ob- 
ject, even  if  this  is  an  active  device  contained  in  some  kind  of  measurement  fixture.  In 
this  way,  all  system  errors  are  removed,  including  those  introduced  by  the  fixture. 

The  position  of  the  reference  plane  is  immaterial,  as  long  as  it  is  kept  constant 
during  calibration  and  measurement.  Also,  the  given  solution  of  the  relevant  network 
equations  is  an  exact,  explicit  solution.  Finally,  the  procedure  is  essentially  independent 
of  the  specific  system  configuration. 
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VI.     CONCLUSION 

Various  sources  have  been  cited  for  the  extraction  of  the  scattering  parameters  of  a 
device  by  removal  of  the  error  and  contributions  of  the  fixturing  setup  [Refs.  8,9,10,11]. 
In  Chapters  II  and  III,  the  de-embedding  procedure  using  microstrip  line  within  the 
given  conditions  was  discussed.  This  project  included  various  aspects  such  as  the 
microstrip  line  fixturing  setup  and  the  construction  and  handling  of  small  experimental 
boards.  But  in  spite  of  the  limitations  of  the  early-model  Network  analyzer  HP  8409B 
which  was  available  for  this  work,  the  best  possible  results  were  obtained  using  repeated 
measurements  on  the  same  experimental  board  to  obtain  mean  values. 

To  extract  the  scattering  parameters  of  the  device  51  ohm  resistor,  the  nine  basic 
parameters  as  discussed  earlier  in  Chapter  IV  (Eqs.  4-11  to  4-19)  were  used.  These 
equations  required  the  scattering  parameters  of  both  fixtures  since  they  are  not  sym- 
metrical. That  is,  individual  measurements  and  calculations  were  required  for  both  fix- 
tures. The  matrix  equation  (Eq.  2-17)  involved  the  three  unknown  scattering  parameters 
Sn,  S22  and  A  =  (SUS22  —  Sl2S2l)  for  fixture  A  and  B.  That  meant  each  fixture  required 
three  equations  for  the  matrix  solution.  Three  different  measurements  and  calculations 
of  microstrip  line  were  used  to  set  up  the  above  matrix  entries,  i.e,  three  different  meas- 
ured reflection  coefficients  from  the  network  analyzer  with  respect  to  the  three  different 
microstrip  loads,  respectively  Tml,  rin2  and  r,n3  and  three  different  calculated  reflection 
coefficients  at  the  terminals  of  the  three  different  loads  rL[,  rL2  and  r\3  were  required. 
The  computed  value  for  Sn  and  S22  are  combined  with  a  knowledge  of  A  which  can  be 
used  to  determine  the  Sl2S2l  product.  Since  Su,  and  S21  for  the  fixtures  always  appear  in 
product  form,  the  magnitude  of  these  terms  can  be  calculated  as  the  square  root  of  the 
product  term  in  the  expression  for  A  . 

The  de-embedded  scattering  parameters  of  the  device  (51  ohm  resistor)  were  found 
through  experiment  to  be  close  to  the  theoretically  calculated  values,  shown  in  magni- 
tude plots  in  Figs.  19,  20,  and  21  and  as  tabulated  in  Table  4.  As  may  be  seen  in  Figs. 
19  to  21,  the  theoretically  calculated  scattering  coefficients  of  the  51 -Ohm  resistor  device 
showed  only  a  slight  continuous  variation  over  the  entire  frequency  range  of  interest. 
The  measured  values  of  the  embedded  device  showed,  however,  rapid  irregular  fluctu- 
ations with  frequency,  reflecting  the  probable  occurrence  of  internal  resonances  within 
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parasitic  discontinuities  of  the  coupling  circuitry.  It  is  the  task  of  the  de-embedding 
procedure  to  remove  these  effects  from  the  measurement.  The  curves  numbered  (3)  in 
the  above  Figures  show  the  de-embedded  measurements,  in  comparison  with  the  the- 
oretical curves.  As  may  be  seen,  the  de-embedded  results  followed  the  expected  values 
closely,  showing  fluctuations  which  are  assumed  to  be  due  to  measurement  error,  and 
network-analvzer  calibration  errors. 


Table  4.     THE  COMPARISON 
OHM  RESISTOR  IN 


OF  THE  OUTPUT  FOR  THE  51 
DB  UNITS 


Freq.(MHz) 

Theoretical 

De- 
embedded 

Embedded 

2000.0 

6.11 

6.35 

11.14 

2300.0 

5.54 

5.12 

10.18 

2600.0 

5.02 

4.22 

7.13 

2900.0 

4.55 

3.87 

6.21 

3200.0 

4.13 

3.19 

7.33 

3500.0 

3.76 

3.16 

11.29 

3800.0 

3.42 

3.91 

11.59 

4100.0 

3.12 

2.66 

6.27 

4400.0 

2.86 

4.11 

6.07 

4700.0 

2.62 

2.25 

9.91 

5000.0 

2.41 

2.11 

11.95 

S12  and  S2l  terms  for  fixture  A  and  B  always  appears  in  product  form,  considering 
their  phase  is  not  important.  As  long  as  the  phase  angles  of  S21  for  fixture  A  and  S21  for 
fixture  B  are  approximately  the  same,  one  can  choose  either  of  the  phase  angles.  Some 
confusion  may  arise,  if  the  scattering  parameters  of  both  fixture  A  and  B  are  specified 
independently  of  the  final  device  (51  ohm  resistor)  parameters.  In  this  case,  the  ap- 
proximate phase  shift  through  one  fixture  can  be  computed,  or  determined  exper- 
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imentally,  and  used  to  choose  the  proper  phase  of  S21  for  fixture  A  and  S21  for  fixture 
B. 

In  Chapter  V,  the  de-embedding  procedure  with  symmetrical  fixtures  was  carried 
out,  using  as  the  device  a  commercial  bandpass  filter  (8-12  GHz).  Standard  hardware 
couplers  were  added  to  represent  the  fixtures.  Since  the  fixtures  were  symmetrical,  the 
de-embedding  procedure  was  relatively  straightforward.  As  may  be  seen  in  Figs.  26  to 
2S.  the  de-embedded  scattering  parameters  of  the  filter  were  very  close  to  the  expectation 
values,  i.e.  of  the  scattering  parameters  of  filter  without  fixtures.  It  may  also  be  noted 
that  the  measured  S-parameters  of  the  embedded  filter  (fixtures  plus  filter)  were  very 
similar  to  those  of  the  filter  alone,  except  of  course  for  phase;  indicative  of  well-matched 
coupling  hardware.  Therefore,  in  this  task,  the  de-embedded  results  were  found  to  be  in 
good  agreement  with  the  expected  values,  as  shown  in  Figs.  26  through  28. 
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APPENDIX  A.     THE  OUTPUTS  OF  THE  NETWORK  ANALYZER  AND 
THEIR  REFLECTION  COEFFICIENTS  FOR  FIXTURE  A  AND  B 


urv's  '  mr- 


Figure  33.       Measurement  of  the  Load  1  (Fixture  A) 


Table  5.       THE  MEASUREMENT  OF  THE  LOAD  HFIXTUF 

l£  A) 

Freq.(MHz) 

(IB  (loss) 

Anele(dec. ) 

r\ 

2000.0 

3.07 

">")  7 

0.65  + jO 

-i- 

2I50.U 

3.S2 

-63.1 

0.29  -  jO 

57 

23<)0.<t 

3.79 

-153.9 

-0.58  -  jo 

28 

2450.0 

3.46 

126.4 

-0.40  4- j() 

54 

2600.0 

3.2o 

60.5 

0.34  -HO 

60 

2750.0 

3 . 3  j 

-1.7 

0.6S  -  jO 

02 

29oo.o 

4.04 

-69.9 

o.22  -  jo 

59 

3050.0 

5.39 

-1-2.8 

-0.4S  -  jo 

25 

3200.O 

6.17 

107.6 

-0.15  +  jo 

4" 

3350.0 

5.32 

17.1 

0.52  + jo 

16 

3500.0 

4.65 

-51.0 

0.37  -  jo 

45 

3650.0 

3.97 

-111.9 

-0.24  -  JO 

59 

3800.0 

4.89 

-175.9 

-0.57  -  jO 

04 

3950.0 

7.73 

104.5 

-0.10  +  jO 

40 

4100.0 

7.36 

-17.1 

0.41  -  jO 

13 

4250.0 

6.70                             -121.4 

-0.24  -  jo 

39 

4400.0 

4.59                             170.4 

-0.58  + jO 

10 

4550.0 

4.12 

114.2 

-0.26 -'-jO 

^7 

4700.0 

4.88 

58.9 

0.29  + jo 

49 

4S50.0 

7.91 

-9.1 

0.40  -  jO 

06 

5000.0 

8.42 

-137.9 

-0.28  -  jO 

25 
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Figure  34.       Measurement  of  the  Load  2  (Fixture  A) 


Table  6.      THE  MEASUREMENT  OF  THE  LOAD  2  (FIXTURE  A> 


Freq.(MHz) 

dB  (loss) 

Ang!e(deg.) 

r.„ 

2O00.0 

2.2S 

9S.2 

-O.II+jO.76 

2150.0 

2.S9 

29.5 

0.62  +  JO.3.5 

230O.O 

3.77 

-49.2 

0.42  -  jO.49 

2450.0 

4.12 

-138.0 

-o.4b  -  ]0.42 

2600.0 

3.75 

136.6 

-0.47^)0.45 

2750.0 

3.31 

65.1 

0.29  +  JO.62 

2900.0 

3.40 

4.9 

0.67-t-ji).06 

3050.0 

3.60 

-54.2 

0.39  -  jO.54 

3200.0 

4.55 

-119.0 

-o.29  -  jO.52 

3350.0 

5.7S 

159.8 

-0.4S-jO.lS 

3500.0 

6.39 

61.8 

U.23  +  J0.42 

3650.0 

5.76 

-25.5 

0.47  -  jO.22 

3800.0 

5.70 

-89.0 

0.01  -  jO.52 

3950.0 

4.69 

-134.2 

-0.41  -  jo.42 

4100.0 

3.91 

173.7 

-0.63  +  jO.07 

4250.0 

4.48 

117.4 

-0.27-jO.53 

4400.0 

6.84 

47.9 

0.31+J0.34 

4550.0 

S.48 

-61.4 

o.lS  -  jo.33 

4700.0 

7.57 

-174.8 

-0.42  -  ]0.04 

4850.0 

4.92 

115.1 

-0.24-«-jU.51 

5000.0 

3.99 

62.0 

0.30  + jo. 56 
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Figure  35. 


Measurement  of  the  Load  3  (Fixture  A) 


Table  7.       THE  MEASUREMENT  OF  THE  LOAD  3  (FIXTURE  A) 

Frcq.(MHz) 

dB  (loss)                    Anele(deg.) 

r,„ 

2noo.o 

2.14                              175.1 

-0.78  +  jO.O? 

2150.0 

2.21 

114.5 

-0.32  +  J0.7I 

2300.0 

2.57 

55.4 

0.42  +  JO.61 

2450.0 

3.25 

-5.1 

0.69  -  jO.06 

2600.0 

3.86 

-74.3 

0.1"  -  J0.62 

2750.0 

4.44 

-156.0 

-0.55  -  jO.24 

2900.0 

4.42 

119.8 

-0.30-r  jO.52 

3050.0 

3. 98 

46.6 

0.43-»-jo.46 

3200.0 

3.89 

-14.2 

0.62  -  jo. in 

3350.0 

3.89 

-66.5 

0.2^  -  jo.59 

3500.0 

4.15 

-120.7 

-0.32  -  jO.53 

3650.0 

5.02 

175.9 

-0.56  +  JO.O4 

3S00.0 

8.22 

9M.0 

-0.06  +  J0.38 

3950.0 

8.95 

9.5 

0.35  +  J0.06 

4100.0 

6.37 

-79.5 

0.09  -  jO.47 

4250.0 

4.35 

-147.2 

-0.51  -  jO.33 

4400.0 

3.49 

157.8 

-0.62-jO.25 

4550.0 

3.37 

108.4 

-0.21  -t-j0.64 

4700.0 

3.90 

58.8 

0.33  +  jo.55 

4850.0 

5.55 

1.4 

0.53  +  jO.OI 

5000.0 

6.59 

-82.6 

0.06  -  jO.46 

59 


JL 
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Figure  36.       Measurement  of  the  Load  1  (Fixture  B) 


Table  8.       THE  MEASUREMENT  OF  THE  LOAD  1  (FIXTURE  B) 

Freq.(MHz) 

dll  (loss) 

Anglei  cleg.) 

r„ 

2000.0 

1.97 

-ir.s 

-0.37  -  ji).71 

2150.0 

1.S9 

-T6.2 

-0.S0  -  jO.05 

2300.0 

2.oo 

13o.2 

-0.51  +JO.01 

2450.0 

2.27 

81 .9 

0.11  i-jO.76 

2600.0 

2.4o 

34.  S 

0.62  +  J0.43 

2750.0 

2.62 

-13.9 

0.72  -  jo. IS 

2900.0 

3.04 

-6". 3 

0.27  -  ju.65 

30.50.0 

3.S7 

-130.3 

-0.41  -  jO.49 

3200. o 

5.05 

151.6 

-0.49  +  JO.27 

3350.0 

5.41 

64.0 

0.24  +  JO.48 

3500.0 

4.91 

-11.7 

0.50  -  jO.  1 2 

3650.0 

3.87 

-67.6 

0.24  .  jo.59 

3SOO.0 

2.90 

-118.9 

-0.35  -  jO.63 

3950.0 

2.55 

-160.3 

-0.72  -  JO. IS 

4100.0 

2.68 

147.0 

-0.62  +  10.40 

4250.0 

3.33 

96.3 

-0.07  +  JO.68 

4400.0 

4.95 

30.4 

0.40  4- ,0.34 

4550.0 

7.33 

-4_v4 

0.31  -  ,0.30 

4700.0 

7.51 

-139.7 

-0.32  -  jo.27 

4850.0 

5.35 

143.0 

-0.434-jO.32 

5000.0 

3.67 

87.9 

0.02  +  JO.65 

6d 
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Figure  37.       Measurement  of  the  Load  2  (Fixture  B) 


Table  9.       THE  MEASUREMENT  OF  THE  LOAD  2  (FIXTURE  R) 

Freq.l  MHz) 

dB  (loss) 

.\ncle(deg.) 

r,„ 

2000.0 

2.32 

-26.9 

0.68  -  jO.35 

2150.0 

2.46 

-91.9 

-O.02  -  jO.75 

2300.0 

2.44 

-155.5 

-0.69  -  jo.31 

2450.0 

2.51 

147.3 

-0.63~jO.4n 

2600.0 

2.41 

96.6 

-O.09-f-jo.75 

2750.0 

2.25 

50. 1 

0.50  + jo. 5y 

2900.0 

2.35 

6.9 

0.76  +  JO.09 

3050.0 

2. 28 

-35.9 

0.62  -  jO.45 

3200.0 

2.51 

-7S.5 

0.15  -  jO.73 

3350.0 

2.69 

-123.9 

-0.41  -  jO.61 

350o.o 

3.25 

-175.1 

-0.69  -  J0.06 

3650.0 

4.27 

126.2 

-0.36-Hjo.4g 

38OO.0 

5.24 

56.1 

0.31  +J0.45 

3950.0 

5.38 

-20.6 

(J. 50  -  jO.19 

4100.0 

4.49 

-90.2 

0.00  -  ]0.6o 

4250.0 

3.49 

-147.9 

-0.57  -  j0.36 

4400.0 

2.76 

163.0 

-o.7o-jo.21 

4550.0 

2.43 

118.9 

-0.37-f-jo.66 

4  700.  <J 

2.36 

77.5 

0.16  +  JO.74 

4850.0 

2.55 

36.2 

0.60  +  J0.44 

5000.0 

3.32 

-8.6 

0.67  -  jO.10 

61 
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Figure  38.       Measurement  of  the  Load  3  (Fixture  B) 


Table   10.      THE  MEASUREMENT  OF  THE  LOAD  3  (FIXTURE  B) 

Freq.(MHz) 

dli  (loss) 

Anele(dec.) 

r 

2000.0 

1.41 

78.5 

0 

I74-J0.83 

2 1 50.0 

1.71                                28.1 

0 

72-J0.39 

23oo.  0 

2.16 

-2o.2 

0 

"0  .  ,0.34 

245o.o 

2.4S 

-84.3 

0 

0-  -  jo.75 

2600.0 

2.69 

-145.9 

-0 

61  -  i0.41 

2750.0 

2.85 

152.4 

-0 

64  ■*-  jo.33 

2900.0 

2.85 

44.8 

-0 

06 +-10.72 

3050.0 

2.75 

43.2 

O 

53  ■+■  jo.So 

3  20O.0 

-»  7- 

-2.s» 

(1 

7">  -  ;0.(l4 

3350.0 

2.08 

-44.3 

0 

53-J0.51 

3500.0 

2.68 

-82.4 

0 

|M    .    ,0.71 

3650.0 

1.93 

-121.1 

-0 

41  -  jO.69 

3S<>0.0 

1.91 

-162.1 

-0 

76  -  10.25 

395o.o 

2.02 

157.0 

-0 

73  + jo.  31 

4  IO0.0 

2.42 

113.2 

-0 

30  +  J0.70 

4250.0 

3.21 

64.h 

0 

3o^j0.62 

4400.0 

4.49 

7.5 

0 

59  +  jo.OS 

4550.0 

5.70 

-60.8 

0 

25  -  ]0.45 

4700.0 

5.8o 

-134. si 

-0 

36  -  jO.36 

4850.0 

4.S9                              158.8 

-O 

53-J0.21 

5000.0 

3.72                             105.1 

-0 

17  +  JU.63 
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APPENDIX  B.     THE  S-PARAMETERS  OF  THE  FIXTURES 


Table   11.      THE  S  -  PARAMETERS  OF  FIXTURE  A 


Freq.(MHz) 

Au 

A 

12 

A 

2 

2000.0 

-0.130  +  J0.0S2 

0.837 

-   jO.143 

0.132 

4-  jO.029 

2150.0 

-0.064  +  j0. 153 

0.668 

-   jO.499 

0.139 

-  jO.019 

2300.0 

0.026  +  j0. 165 

0.355 

-   jO.738 

0.130 

-  jO.OSO 

2450.0 

0.099  +  J0.125 

0.021 

+  jO.802 

0.100 

-  jO.  147 

2600.0 

0.141  +  J0.069 

0.380 

+  jO.701 

0.034 

-  jO.178 

2750.0 

0.151    -  jO.007 

0.653 

+■  jO.443 

-0.018 

-  jO.180 

2900.0 

0.091    -  ]0.063 

0.777 

+  jO.087 

-0.032 

-  jO.128 

3050.0 

0.048   -   jO.055 

0.717 

-  J0.287 

-0.003 

-  jO.074 

3200.0 

0.042   -  J0.026 

0.485 

-   jO.576 

0.060 

-  jO.037 

3350.0 

0.098   -   jO.025 

0.151 

-  jO.720 

0.180 

-  jO.031 

3500.0 

0.137   -  J0.113 

0.199 

+  jO.679 

0.2S8 

-  jO.088 

3650.0 

0.054   -  jO.265 

0.470 

+  jO.482 

0.330 

-  jO.184 

3S00.0 

-0.148   -  J0.261 

0.579 

+  jO.214 

0.327 

-   jO. 254 

3950.0 

-0.258   -  jO.163 

0.625 

-  jO.036 

0.321 

-  J0.239 

4100.0 

-0.334  +  J0.061 

0.581 

-  jO.305 

0.299 

-  jO.322 

4250.0 

-0.256  +  jO.288 

0.318 

-  jO.506 

0.442 

-   jO.181 

4400.0 

0.174  +  J0.317 

0.008 

-   jO.477 

0.713 

-   jO.044 

4550.0 

-0.502   -  jO.872 

0.225 

-   jl.580 

1.743 

-   jl.374 

4700.0 

-0.064   -  jO.531 

0.670 

-  jO.309 

0.034 

-   jO.761 

4850.0 

0.183  +  J0.229 

0.651 

+■  jO.057 

0.137 

-   jO.304 

5000.0 

0.062  +  jO.136 

0.646 

-  jO.289 

0.219 

-   jO.013 
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Table  12.      THE  S  -  PARAMETERS  OF  FIXTURE  B 


Frcq.(MHz) 

B» 

B 

2 

B22 

2000.0 

-0.060  4-  jO.  107 

0.342 

-  jO.834 

0.036   -  jO.076 

2150.0 

-0.007  4-  J0.116 

0.057 

-   J0.S92 

0.008   -   J0.081 

2300.0 

0.027  +  jO.  100 

0.230 

4-  jO.850 

-0.032   -  jO.085 

2450.0 

0.049  4-  J0.091 

0.475 

+  jO.723 

-0.083   -  j.082 

2600.0 

0.0S0  4-  jO.084 

0.669 

+  jO.533 

-0.115   -  jO.075 

2750.0 

0.124  4-  jO.063 

0.792 

-  jO.291 

-0.125   -  jO.062 

2900.0 

0.152  4-  jO.007 

0.831 

4-  jO.028 

-0.136   -  jO.029 

3050.0 

0.158   -  jO.073 

0.788 

-  jO.229 

-0.125   -  jO.007 

3200.0 

0.106   -   jO.143 

0.669 

-  J0.451 

-0.110  4-  jO.016 

3350.0 

0.026   -   jO.168 

0.494 

-  jO.630 

-0.064  4-  J0.036 

3500.0 

-0.035   -  jO.149 

0.268 

-   jO.741 

0.016  4-  J0.059 

3650.0 

-0.092   -  jO.  112 

0.033 

-  jO.790 

0.108  4-  J0.019 

3800.0 

-0.076  -  jO.072 

0.230 

4-  J0.762 

0.165  4-  J0.018 

3950.0 

-0.069   -  jO.097 

0.468 

4-  J0.633 

0.261  4-  J0.019 

4100.0 

-0.137   -  J0.157 

0.636 

+  J0.412 

0.389  4-  J0.014 

4250.0 

-0.331    -  J0.130 

0.681 

4-  jO.135 

0.543  4-  jO.018 

4400.0 

-0.572  +  jO.161 

0.478 

-  jO.116 

0.808  4-  jO.024 

4550.0 

-0.114  4-  jl.171 

1.062 

4-  J0.589 

1.625   -  jO.527 

4700.0 

0.230  -  jO.171 

0.577 

-  jO.759 

0.163   -  jO.968 

4850.0 

-0.065  4-  jO.166 

0.306 

-  J0.676 

0.219   -  jO.546 

5000.0 

0.124  4-  J0.326 

0.091 

-  jO.685 

0.291    -  jO.446 
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APPENDIX  C.     THE  S-PARAMETERS  OF  THE  EMBEDDED  AND 
DE-EMBEDDED  51  OHM  RESISTOR 


Table  13.     THE  THEORETICAL  S  -  PARAMETERS  OF  THE  51 
OHM  RESISTOR  IN  DB  UMTS. 


Freq.(MHz) 

Su 

su 

S2] 

si2 

2000.0 

6.11 

4.27 

4.27 

6.11 

2150.0 

5.82 

4.37 

4.37 

5.82 

2300.0 

5.54 

4.47 

4.47 

5.54 

2450.0 

5.28 

4.58 

4.58 

5.28 

2600.0 

5.02 

4.69 

4.69 

5.02 

2750.0 

4.78 

4.81 

4.81 

4.78 

2900.0 

4.55 

4.93 

4.93 

4.55 

3050.0 

4.34 

5.05 

5.05 

4.34 

3200.0 

4.13 

5.17 

5.17 

4.13 

3350.0 

3.94 

5.30 

5.30 

3.94 

3500.0 

3.76 

5.43 

5.43 

3.76 

3650.0 

3.58 

5.56 

5.56 

3.58 

3S00.0 

3.42 

5.69 

5.69 

3.42 

3950.0 

3.27 

5.82 

5.82 

3.27 

4100.0 

3.12 

5.95 

5.95 

3.12 

4250.0 

2.99 

6.09 

6.09 

2.99 

4400.0 

2.86 

6.22 

6.22 

2.86 

4550.0 

2.73 

6.36 

6.36 

2.73 

4700.0 

2.62 

6.49 

6.49 

2.62 

4S50.0 

2.51 

6.63 

6.63 

2.51 

5000.0 

2.41 

6.76 

6.76 

2.41 
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Table  14.      THE  S  -  PARAMETERS  OF  BOTH  FIXTURE  AND  51 
OHM  RESISTOR 


Freq.(MHz) 

Mu 

Mn 

A/21 

Ma 

2000.0 

11.14 

4.40 

4.60 

6.45 

2150.0 

11.38 

4.64 

4.71 

6.27 

2300.0 

10.18 

4.90 

4.80 

6.05 

2450.0 

8.47 

5.36 

5.35 

5.86 

2600.0 

7.13 

5.87 

5.84 

5.51 

2750.0 

6.35 

6.30 

6.25 

4.99 

2900.0 

6.21 

6.99 

6.99 

4.68 

3050.0 

6.51 

7.54 

7.49 

4.48 

3200.0 

7.33 

7.32 

7.25 

4.43 

3350.0 

8.92 

7.04 

6.65 

4.94 

3500.0 

11.29 

7.10 

5.81 

5.75 

3650.0 

12.01 

6.46 

6.31 

6.18 

3800.0 

11.59 

7.30 

7.25 

6.91 

3950.0 

9.15 

6.45 

6.90 

7.83 

4100.0 

6.27 

6.35 

6.56 

9.07 

4250.0 

5.53 

6.00 

6.25 

1.78 

4400.0 

6.07 

5.61 

5.85 

12.72 

4550.0 

7.54 

5.30 

5.50 

12.53 

4700.0 

9.91 

5.34 

5.39 

10.00 

4850.0 

10.46 

5.70 

5.66 

7.82 

5000.0 

11.95 

6.10 

6.05 

6.53 
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Table  15.       THE  DE  -  EMBEDDED  OUTPUT  OF  51  OHM  RE- 
SISTOR DEVICE  UNDER  TEST. 


Freq.(MHz) 

Su 

sl2 

s2i 

Sr_ 

2000.0 

6.35 

3.62 

3.82 

4.92 

2150.0 

5.72 

3.97 

3.64 

5.84 

2300.0 

5.12 

3.97 

3.97 

4.74 

2450.0 

4.57 

3.18 

3.17 

3.75 

2600.0 

4.22 

3.39 

3.36 

3.04 

2750.0 

3.61 

3.66 

3.61 

4.74 

2900.0 

3.87 

4.27 

4.27 

3.77 

3050.0 

3.64 

4.91 

4.86 

2.91 

3200.0 

3.19 

4.51 

4.44 

2.67 

3350.0 

3.03 

4.16 

3.77 

2.41 

3500.0 

3.16 

4.72 

4.43 

3.79 

3650.0 

3.62 

4.65 

4.50 

4.45 

3800.0 

3.91 

4.56 

4.51 

3.72 

3950.0 

3.35 

5.38 

4.83 

2.73 

4100.0 

2.66 

5.06 

5.27 

4.66 

4250.0 

3.45 

5.71 

5.46 

5.31 

4400.0 

4.11 

5.83 

6.07 

2.52 

4550.0 

3.17 

6.36 

6.56 

2.07 

4700.0 

2.25 

6.73 

6.68 

3.34 

4850.0 

2.08 

6.73 

6.69 

2.95 

5000.0 

2.11 

6.96 

6.91 

2.02 
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APPEiNDIX  D.     THE  S-PARAMETERS  OF  THE  FIXTURE,  FILTER 
EMBEDDED,  FILTER  WITHOUT  FIXTURES  AND  FILTER 

DE-EMBEDDED 


Table   16.      THE  S 

-  PARAMETERS  OF  THE  FIXTURE 

Froq.(MHz) 

s„ 

su 

sa 

6000.0 

-0.08  4-  jO.17 

0.81  +  jO.55 

-0.09  +  J0.13 

6320.0 

-0.09  4-  jO.04 

0.76   -   jO.61 

-0.17  4-  jO.10 

6640.0 

-0.06  +  jO.04 

-0.39   -   jO.90 

-0.03  4-  J0.12 

6960.0 

-0.05  +  J0.03 

-0.98  4-  jO.16 

-0.08  4-  J0.07 

72S0.0 

-0.19  +  J0.03 

-0.08  4-  jO.97 

-0.18  4-  jo.15 

7600.0 

-0.13  +  jO.14 

0.92  4-  jO.36 

-0.05  4-  jO.18 

7920.0 

0.01  4-  jO.07 

0.57   -  jO.81 

0.01  +  j0.ll 

8240.0 

-0.07  +  jO.OO 

-0.63   -  jO.76 

0.02  4-  J0.11 

S560.0 

-0.12  +  jO.04 

-0.90  4-  jO.42 

0.03   -  jO.01 

8880.0 

-0.10  4-  J0.07 

0.18  4-  J0.97 

-0.11    -  jO.03 

9200.0 

0.00  +  jo.19 

0.96  4-  J0.06 

-0.03  4-  J0.20 

9520.0 

0.10  4-  jO.14 

0.33   -   jO.91 

0.11  4-  J0.11 

9S40.0 

0.11    -  jO.04 

-0.81    -   jO.55 

-0.02   -   J0.09 

10160.0 

0.12   -  jO.19 

-0.74  4-  jO.66 

-0.03    -   J0.03 

10480.0 

-0.03   -  J0.13 

0.42  4-  J0.88 

0.03  4-  J0.02 

10800.0 

-0.02   -  jO.14 

0.96   -  jO.18 

-0.02   -  J0.09 

11120.0 

-0.18   -  jO.12 

0.06   -  J0.98 

-0.16  +  J0.03 

11440.0 

-0.14   -  J0.03 

-0.91    -  jO.33 

-0.08  4-  J0.02 

11760.0 

-0.02   -  jO.14 

-0.57  +  J0.79 

-0.04   -  jO.04 

12080.0 

-0.04   -  jO.04 

0.63  +  J0.75 

0.04  +  J0.06 

12400.0 

-0.02   -  jO.10 

0.S6  -  ]0.45 

-0.02   -  jO.04 
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Table  17.      THE  S 

-  PARAMETERS  OF  FILTER  +  FIXTURE 

Freq.(MHz) 

s„ 

sa 

s:2 

6000.0 

-0.79   -  J0.45 

0.00  +  jO.OO 

0.54  4-  jO.75 

6320.0 

0.S9   -   J0.14 

0.00  4-  jO.OO 

-0.36  -  J0.S2 

6640.0 

-0.65  +  jO.62 

0.01  4-  jO.01 

-0.02  4-  jO.90 

6960.0 

0.11    -   jO.91 

0.01  +  J0.01 

0.05   -  jO.89 

72S0.0 

0.62  4-  jO.64 

0.01  4-  jO.01 

-0.24  4-  J0.90 

7600.0 

-0.S0  4-  J0.45 

-0.02  +  jO.02 

0.02   -   jO.91 

7920.0 

0.11    -   J0.24 

-0.30  4-  J0.86 

0.00  4-  J0.25 

8240.0 

-0.13  4-  jO.20 

-0.34  +  jO.86 

-0.01    -  jO.15 

8560.0 

-0.01  4-  jO.36 

-0.82  4-  jO.24 

0.15  4-  J0.27 

SSSO.O 

0.01  4-  J0.04 

-0.61    -  J0.69 

-0.01    -   jO.05 

9200.0 

0.31  +  jO.01 

0.56   -  jO.65 

0.03  4-  J0.31 

9520.0 

0.27  4-  J0.29 

0.76  4-  jO.41 

-0.28   -  J0.06 

9840.0 

-0.16  4-  J0.22 

-0.48  4-  J0.76 

0.07   -  jO.23 

10160.0 

-0.06  +  jO.23 

-0.73   -  jO.50 

0.00  4-  jO.29 

104S0.0 

-0.06   -  jO.26 

0.47   -  J0.76 

-0.26   -  jO.03 

10800.0 

0.21    -  jO.02 

0.79  4-  J0.46 

-0.05   -  jO.19 

11120.0 

-0.17   -  J0.07 

-0.45  4-  J0.78 

-0.14  -  jO.02 

11440.0 

-0.14  +  jO.OO 

-0.83   -  jO.37 

0.14  4-  jO.01 

11760.0 

0.19  +  jO.07 

0.16   -  jO.86 

0.17  4-  jO.05 

12080.0 

0.25   -   jO.24 

0.72   -  jO.28 

-0.48   -  J0.13 

12400.0 

-0.15  4-  J0.70 

0.74  4-  jO.18 

0.59  4-  J0.12 
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Table  18.      THE  S  -  PARAMETERS  OF  FILTER  WITHOUT  FIX- 
TURES IN  DB  UNITS 


Freq.(MHz) 

S» 

si2 

s2i 

s22 

6000.0 

1.57 

83.91 

88.79 

0.62 

6320.0 

1.96 

91.38 

105.37 

1.69 

6640.0 

0.47 

75.81 

78.12 

0.19 

6960.0 

0.87 

69.94 

69.39 

0.85 

72S0.0 

0.77 

54.40 

51.94 

0.86 

7600.0 

0.70 

29.99 

31.60 

1.12 

7920.0 

10.17 

0.57 

0.50 

9.15 

S240.0 

11.32 

0.15 

0.16 

12.00 

8560.0 

8.62 

0.45 

0.45 

10.05 

SS80.0 

15.77 

0.20 

0.08 

13.81 

9200.0 

9.66 

0.40 

0.20 

9.88 

9520.0 

10.97 

0.25 

0.33 

11.81 

9840.0 

9.65 

0.27 

0.66 

9.11 

10160.0 

8.39 

0.45 

0.40 

10.81 

10480.0 

12.78 

0.15 

0.20 

12.01 

10S00.0 

11.35 

0.43 

0.45 

11.87 

11120.0 

13.64 

0.30 

0.20 

14.39 

11440.0 

10.28 

0.31 

0.46 

16.09 

11760.0 

15.44 

0.60 

0.41 

15.84 

12080.0 

7.35 

1.40 

1.10 

9.98 

12400.0 

8.85 

1.45 

0.96 

11.51 
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Table  19.      THE  DE-EMBEDDED  OUTPUT  OF  FILTER  UNDER 
TEST 


Freq.(MHz) 

sn 

sa 

s2i 

S22 

6000.0 

1.58 

37.03 

90.64 

0.28 

6320.0 

2.03 

41.18 

91.12 

1.15 

6640.0 

1.11 

41.62 

76.82 

0.14 

6960.0 

2.32 

59.94 

71.12 

0.90 

72S0.0 

0.68 

51.99 

53.34 

1.21 

7600.0 

0.21 

30.06 

31.45 

0.92 

7920.0 

10.50 

1.27 

1.51 

9.64 

S240.0 

10.99 

0.59 

0.60 

13.25 

8560.0 

7.97 

1.25 

1.18 

10.17 

SSSO.O 

12. S3 

1.12 

0.55 

10.50 

9200.0 

9.57 

1.36 

0.70 

10.77 

9520.0 

8.80 

0.35 

0.94 

9.70 

9S40.0 

10.83 

0.11 

0.81 

8.38 

10160.0 

7.58 

0.74 

0.90 

10.05 

104S0.0 

11.96 

0.62 

0.75 

9.76 

10S00.0 

12.39 

0.67 

0.73 

9.99 

11120.0 

11.49 

0.14 

0.20 

9.78 

11440.0 

8.94 

0.38 

0.76 

15.48 

11760.0 

13.73 

0.76 

0.94 

14.58 

12080.0 

8.20 

2.11 

1.41 

6.S5 

12400.0 

8.41 

1.77 

1.38 

7.79 
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APPENDIX  E.     FORTRAN  PROGRAM  FOR  DE-EMBEDDING  THE  51 

OHM  RESISTOR 


DIMENSION  X(200),MDB11(200),PSE11(200) 
MDB12(200),PSE12(200) 
MDB21(200),PSE21(200) 
MDB22(200) ,PSE22(200) ,MDB222(200) 
REAL  MDBU,MDB12,MDB21,MDB22,MDB222 
MAG 1 1 , MAG 1 2 , MAG2 1 , MAG2  2 
IM11,IM12,IM21,IM22 
COMPLEX  S11,S12,S21,S22,A12A21,B12B21 

A11,A12,A21,A22,AA121,A12B12,A21B21 
B11,B12,B21,B22,BB121 
M11,M12,M21,M22 
E,FF,G,H,N,AABB 
PI=4.*ATAN(1.  ) 

DO  10  1=1,21 

F=18.  5*10.  **8+l.  5*10.  **8*FL0AT(I) 
X(I)=F/1000000. 

CALL  FIX1(A11,A12,A21,A22,A12A21,AMAG12,PSEA12,F) 

CALL  TMFS(M11,M12,M21,M22,F) 

CALL  FIX2(B11,B12,B21,B22,B12B21,BMAG12,PSEB12,F) 


TMAG=AMAG12*BMAG12 

TPSE=PSEA12+PSEB12 

XRE=TMAG*COS(TPSE) 

XIM=TMAG*SIN(TPSE) 

A12B12=CMPLX(XRE,XIM) 

A21B21=A12B12 

E=M12/(A12B12) 

FF=M21/(A21B21) 

G=(M11-A11)/A12A21 

H=(M22-B22)/B12B21 

N=( 1.  +A22*G)*( 1.  +B11*H)-A22*B11*E*FF 


S11=(G*(1.+B11*H)-B11*E*FF)/N 

S12=E/N 

S21=FF/N 

S22=(H»(  1.  +A22*G)-A22*E*FF)/N 


QIM1=AIMAG(S11) 
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QRE1=REAL(S11) 

QIM2=AIMAG(S12) 

QRE2=REAL(S12) 

QIM3=AIMAG(S21) 

QRE3=REAL(S21) 

QIM4=AIMAG(S22) 

QRE4=REAL(S22) 

FIND  A  MAGNITUDE  AND  PHASE( ANGLE) 

MAG11=CABS(S11) 

MAG12=CABS(S12) 

MAG21=CABS(S21) 

MAG22=CABS(S22) 

MDB11(I)=ABS(20*ALOG10(MAG11)) 

MDB12(I)=ABS(20*ALOG10(MAG12)) 

MDB21(I)=ABS(20*ALOG10(MAG21)) 

MDB22(I)=ABS(20*ALOG10(MAG22)) 

PSE11(I)=ATAN2(QIM1,QRE1)*180./PI 

PSE12(I)=ATAN2(QIM2,QRE2)*180.  /PI 

PSE21(I)=ATAN2(QIM3,QRE3)*180.  /PI 

PSE22(I)=ATAN2(QIM4,QRE4)*180.  /PI 

WRITE(6,101)  X(I),MDB11(I),MDB12(I),MDB21(I),MDB22(I), 
1  PSE11(I),PSE12(I),PSE21(I),PSE22(I) 

101   FORMAT(F6. 0,8(1X,F7. 2)) 

10    CONTINUE 
END 


Q  ^VVr^VVr•5V'5V'5V')V''^V•)V■5V^VVf■)V'^V'^V')V'!V•)V■lV')V'>V■)V')V•>V')V■)V'!VVf'^V^V**•^V'A, 

SUBROUTINE  TMFS  (Dll ,D12,D21,D22 ,F) 

Q  VrVfVcVc,')VVc,5'r'ilf-<V';'f'>ViV'5V')V'jV'5V,)V'}V'5ViV^!r*V'5V,)V')V,>V'iV')V'5VVf,>V,)V,A'')V">V 

COMPLEX  D11,D12,D21,D22 
PI=4.*ATAN(1.  ) 

INPUT  DATA  FILE  IS  DISK  9 
(  BOTH  FIXTURE  +  DEVICE  ) 

READ(9,21)DB11,ANG11,DB21,ANG21,DB12,ANG12,DB22,ANG22 
21    FORMAT(8(F7.2)) 

Vll=10.**(-DBll/20.  ) 
V12=10.**(-DB12/20.  ) 
V21=10.**(-DB21/20.  ) 
V22=10.**(-DB22/20. ) 
RAD11=PI/180.*ANG11 
RAD12=PI/180.*ANG12 
RAD21=PI/180.*ANG21 
RAD22=PI/180.*ANG22 

QRE11=V11*C0S( RADII) 
QRE12=V12*COS(RAD12) 
QRE21=V21*COS(RAD21) 
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QRE22=V22*C0S(RAD22) 
QIM11=V11*SIN(RAD11) 
QIM12=V12*SIN(RAD12) 
QIM21=V21*SIN(RAD21) 
QIM22=V22*SIN(RAD22) 

D11=CMPLX(QRE11,QIM11) 
D12=CMPLX(QRE12,QIM12) 
D2 1=CMPLX( QRE2 1 , QIM2 1 ) 
D22=CMPLX(QRE22,QIM22) 

RETURN 
END 


SUBROUTINE   FIX1(S11 ,S12 ,S21 ,S22,SS ,SMAG12,SPSE12 ,F) 

FOR  FIXTURE  A 

COMPLEX  A(3,3),H(3,3),B(1,3),C(3),D(3,3),SS121 
COMPLEX  311,512,321,822,33 

CALL  DATA1(A,B,F) 

CALL  CPX(A) 

DO   10    1=1,3 

C(I)=A(1,I)*B(1,1)+A(2,I)*B(1,2)+A(3,I)*B(1,3) 

10  CONTINUE 

S11=C(1) 
S22=C(2) 
SS121=C(3) 
SS=S11*S22-SS121 

SIM=AIMAG(SS) 

SRE=REAL(SS) 

SMAGSS=CABS(SS) 

PSESS=ATAN2(SIM,SRE) 

SPSE12=0.5*PSESS 

SPSE21=SPSE12 

SMAG12=SMAGSS**0.  5 

SMAG21=SMAG12 

S12=(S11*S22-SS121)**0. 5 

S21=S12 

RETURN 

END 


SUBROUTINE  DATA1(D,B,F) 

p  y -  Jm  «' -  tfm  JU  JL  Jm  JU  J*  «J-  -J  -  J*  JU  -'-  J-  Sfc  J-  J-  J-  ^*  - '.-  »'-  V* 

COMPLEX  D(3,3),B(1,3),RIN1,RIN2,RIN3,RL1,RL2,RL3 
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CALL  RL0AD(RL1,RL2,RL3,F) 
CALL  RIN(RIN1,RIN2,RIN3) 


D(l 

,D  = 

=  1. /RL1 

D(l 

,2)  = 

=1.  /RL2 

D(l: 

,3)  = 

=  1.  /RL3 

D(2. 

,1): 

=RIN1 

D(2, 

,2)  = 

=RIN2 

D(2; 

,3)  = 

=RIN3 

D(3: 

,D  = 

=(-i.  ,-0) 

D(3: 

,2)  = 

=(-i.  ,.0) 

D(3, 

,3)  = 

=(-i.  ,.0) 

B(l, 

,1)  = 

=RIN1/RL1 

B(l, 

,2)  = 

=RIN2/RL2 

B(l, 

3)  = 

=RIN3/RL3 

RETURN 
END 

SUBROUTINE  RIN(RIN1,RIN2 ,RIN3) 

COMPLEX  RIN1,RIN2,RIN3 

CALL  F1DAT(DB111,ANG111,DB112,ANG112,DB113,ANG113) 

PI=4.*ATAN(1.  ) 
INPUT  DATA  FILE  IS  DISK8 
(FIXTURE  1  ONLY) 

Vlll=10.**(-DBlll/20.  ) 

V112=10.**(-DB112/20.  ) 

V113=10.**(-DB113/20.  ) 

RAD111=PI/180.*ANG111 

RAD112=PI/180.*ANG112 

RAD113=PI/180.*ANG113 

QRE111=V111*C0S(RAD111) 
QRE112=V112*C0S(RAD112) 
QRE113=V113*C0S(RAD113) 
QIM111=V111*SIN(RAD111) 
QIM112=V112*SIN(RAD112) 
QIM113=V113*SIN(RAD113) 

RIN1=CMPLX(QRE111,QIM111) 
RIN2=CMPLX(QRE112,QIM112) 
RIN3=CMPLX(QRE113,QIM113) 

RETURN 
END 


SUBROUTINE  F1DAT(DB111 , ANG1 1 1 , DB1 12 , ANG112 ,DB113 , ANG113) 
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INPUT  DATA  FILE   IS  DISK  8 
(For  FIXTURE  A) 

READ(8,20)DB111,ANG111,DB112,ANG112,DB113,ANG113 
20  F0RMAT(6(F7.2)) 

RETURN 
END 


r*  y-  y-  y-  y-  y.  y -  y-  yf  y-  y  -  y-  y-  y  -  y  r  y-  y-  y *  y-  y-  -  'ff  y-  y-  y  -  y-  yr  y  -  y-  y-  y-  y  -  y-  y -  y  -  yf  y  -  y  -  y-  y-  y  -  y-  y  -  y  -  y  -  y-  y  r  y-  yf  yf 

SUBROUTINE   FIX2(S11 ,S12,S21 ,S22 ,SSS ,SMAG12,SPSE12 ,F) 

p  J*  -J-  y .  y  -  y .  jl  -.*-  ju  * •  -  y  -  y-  y-  y .  y -  y  .  y-  y  -  «j- «;-  y -  y- y-  y .  y -  y .  y -  y  -  y-  y -  J*  JU  y  -  >J>  y -  y -  - '  -  y -  y -  y.  y-  y-  y  .  y - y-  y  -  y  - y*  y v- 

FOR   FIXTURE    2 

COMPLEX  A(3,3),H(3,3),B(1,3),C(3)JD(3,3),SS121 
COMPLEX   S11,S12,S21,S22,SSS 

CALL  DATA2(A,B,F) 

CALL  CPX(A) 

DO   10    1=1,3 

C(I)=A(l,I)*B(l,l)+A(2,I)*B(l,2)+A(3,I)*B(;i,3) 
10  CONTINUE 

S11=C(1) 
S22=C(2) 
SS121=C(3) 
SSS=S11*S22-SS121 

SIM=AIMAG(SSS) 

SRE=REAL(SSS) 

SMAGSS=CABS(SSS) 

PSESS=ATAN2(SIM,SRE) 

SPSE12=0.5*PSESS 

SPSE21=SPSE12 

SMAG12=SMAGSS**0. 5 

SMAG21=SMAG12 

S12=(S11*S22-SS121)**0.  5 

S21=S12 

RETURN 
END 


SUBROUTINE  DATA2(D,B,F) 

p  y-y-y-yfyry-y-y-y-y-y-y-y-y-y.y-y-y-y-y-y-y-y- 

COMPLEX  D( 3 , 3) ,B( 1 , 3) ,RIN1 ,RIN2 ,RIN3 ,RL1 ,RL2 ,RL3 
CALL  RL0AD(RL1,RL2,RL3,F) 
CALL  R2IN(RIN1,RIN2,RIN3) 


76 


D(1,1)  =  1./RL1 

D(l,2)  =  l.  /RL2 

D(l,3)  =  l.  /RL3 

D(2,1)=RIN1 

D(2,2)=RIN2 

D(2,3)=RIN3 

D(3,l)=(-1.  ,.0) 

D(3,2)=(-l.  ,.0) 

D(3,3)=(-l.  ,.0) 

B(1,1)=RIN1/RL1 

B(1,2)=RIN2/RL2 

B(1,3)=RIN3/RL3 

RETURN 

END 


SUBROUTINE  R2IN(RIN1,RIN2,RIN3) 

ic  iV  it  -i*s  Vr  Vr  Vr  sV  Vr  J-  Vr  *  V  -*-  *■'-  -'-  J-  J-  -'-  J*  J-  ■J-  -'-  *•-  •J-  J-  *  V  J- J-  J-  J-  J- 


c 


COMPLEX  RIN1,RIN2,RIN3 

CALL  F2DAT(DB111,ANG111,DB112,ANG112,DB113,ANG113) 

PI=4.  *ATAN(  1.  ) 

Vlll=10.**(-DBlll/20.  ) 
V112=10.**(-DB112/20. ) 
V113=10.**(-DB113/20.  ) 
RAD111=PI/180.*ANG111 
RAD112=PI/180.*ANG112 
RAD113=PI/180.*ANG113 

QRE111=V111*C0S(RAD111) 
QRE112=V112*COS(RAD112) 
QRE113=V113*COS(RAD113) 
QIM111=V111*SIN(RAD111) 
QIM112=V112*SIN(RAD112) 
QIM113=V113*SIN(RAD113) 

RIN1=CMPLX(QRE111,QIM111) 
RIN2=CMPLX(QRE112,QIM112) 
RIN3=CMPLX(QRE113,QIM113) 

RETURN 
END 


/-i  j-  y-  j-j-jf  jpyf  j.y.yfj*j- y*j^j.y.y.y.  y.y,y-y-y^ 

SUBROUTINE  F2DAT(DB111 , ANG1 1 1 .DB112 ,ANG112,DB113 ,ANG113) 

p         j-  -  ■  -  j .  y . . '-  y ..  u-  j-  y-  jl  j  _  j-  y  -  jj .  y-  y  -  y_  y.  y .  y  f  »'•»'-  -'-  •' -  -';  y-  »'  -  -/-  y-  y-  -'-  -'-  y-  J-  y-  -'-  y-  y-  y-  Vr  *  V  V*  V'  V*  Vr  -,V  V  r  V*  Vc  Vr  V*  ~'r  *  V  ~'r  V  r  v  Vc 

INPUT  DATA  FILE   IS  DISK   10 
(For  FIXTURE   B) 

READ(10,20)DB111,ANG111,DB112)ANG112,DB113,ANG113 
20  FORMAT(6(F7.2)) 
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RETURN 
END 


SUBROUTINE  CPX(A) 

Q  V.-VfiViVVfiVVfVfjV'jV'iV'jV'sVy.-iV'j'oV 

INTEGER   IP( 3) 

COMPLEX  A(3,3),W(3),DET(2) 

CALL  CGEFA(A,3,3,IP,INF0) 

CALL  CGEDI(A,3,3,IP,DET,W,01) 

RETURN 
END 


ri         y.  y-  y-  y-  y-  y-  y-  y-  j-  yf  ye  y*  y*  y*  y  *  y .  y-  y-  y*  yf  yf  y-  y.  yf  y*  yf  jl  yf  y-  y-  yp 

SUBROUTINE  RL0AD(RL1 ,RL2 ,RL3 ,F) 

p  y  -  y  .  y-  -i-  y-  y  .  y  -  y -  y -  y  -  y-  y-  y  -  y-  y  -  y-  y *  y-  y  -  -'-  y  -  y .  y  .  •£.  ju  ju  ju  y,  j1.  jl  -j. 


COMPLEX  RL1 , RL2 , RL3 , ZO , ZL1 , ZL2 , ZL3 

REAL  L1,L2,L3 

PI=4.  *ATAN(  1.  ) 

C=3.*10**8 

Z0S=68.  2 

ERE=2.  833 

ERES=(ERE)**0.5 

L1=0. 04 

L2=0.  03 

L3=0.  02 

A1=C0S(2.  *PI*F*L1*ERES/C) 

B1=SIN(2.*PI*F*L1*ERES/C) 

A2=C0S(2. *PI*F*L2*ERES/C) 

B2=SIN(2. *PI*F*L2*ERES/C) 

A3=C0S(2. *PI*F*L3*ERES/C) 

B3=SIN(2.*PI*F*L3*ERES/C) 

R1=-Z0S*(A1/B1) 

R2=-Z0S*(A2/B2) 

R3=-Z0S*(A3/B3) 

ZL1=CMPLX(0. ,R1) 

ZL2=CMPLX(0. ,R2) 

ZL3=CMPLX(0. ,R3) 

ZO=CMPLX(50.  ,0.  ) 

RL1=(ZL1-Z0)/(ZL1+Z0) 

RL2=(ZL2-Z0)/(ZL2+Z0) 

RL3=(ZL3-Z0)/(ZL3+Z0) 

RETURN 

END 


y m  »v  •.'-  -J-  --'  -  -/.»  «■*—  •.'*  y-  «j-  «J--  — '-  -.'-  >x>  y*  —'•  .j»  «.'-  y*  «j *  »t*  «/-  <j«  .,',»  y^  -j^  -.*..  y-  y^  «.'»  -j« 
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THE  INPUT  DATA  FILE  FOR  DISK  8 
3  DIFFERENT  LOAD  FOR  FIXTURE  A 

-.V-VVf-sV-iV'iVVf-jVVrVcVr-V'jVVfVo'f-iVVf-.VVr-VVfVr-.VVrVfVrV.-VrVfjV 


3.07 

22. 

7 

2. 

28 

98. 

2 

2.  14 

175. 

1 

3.82 

-63. 

1 

2. 

89 

29. 

5 

2.  21 

114. 

5 

3.  79 

-153. 

9 

3. 

77 

-49. 

2 

2.57 

55. 

4 

3.46 

126. 

4 

4. 

12 

-138. 

0 

3.  25 

-5. 

1 

3.  20 

60. 

5 

3. 

75 

136. 

6 

3.  86 

-74. 

3 

3.  33 

-1. 

7 

3. 

31 

65. 

1 

4.44 

-156. 

0 

4.  04 

-69. 

9 

3. 

40 

4. 

9 

4.42 

119. 

3 

5.  39 

-152. 

8 

3. 

60 

-54. 

2 

3.  98 

46. 

6 

6.  17 

107. 

6 

4. 

55 

-119. 

0 

3.  89 

-14. 

2 

5.  32 

17. 

1 

5. 

78 

159. 

8 

3.  89 

-66. 

5 

4.  65 

-51. 

0 

6. 

39 

61. 

8 

4.  15 

-120. 

7 

3.  97 

-111. 

9 

5. 

76 

-25. 

5 

5.  02 

175. 

9 

4.89 

-175. 

9 

5. 

70 

-89. 

0 

8.  22 

99. 

0 

7.  73 

104. 

5 

4. 

69 

-134. 

2 

8.  95 

9. 

5 

7.  36 

-17. 

1 

3. 

91 

173. 

7 

6.  37 

-79. 

5 

6.  70 

-121. 

4 

4. 

48 

117. 

4 

4.  35 

-147. 

2 

4.59 

170. 

4 

6. 

84 

47. 

9 

3.49 

157. 

8 

4.  12 

114. 

2 

8. 

48 

-61. 

4 

3.  37 

108. 

4 

4.  88 

58. 

9 

7. 

57 

-174. 

8 

3.  90 

58. 

S 

7.  91 

-9. 

1 

4. 

92 

115. 

1 

5.55 

1. 

4 

8.42 

-137. 

9 

3. 

99 

62. 

0 

6.59 

-82. 

6 

THE  INPUT  DATA  FILE  FOR  DISK  9 
BOTH  FIXTURE  AND  51  OHM  RES. 


11. 

14 

-9.  3 

4.60 

51.0 

4.40    51.4 

6.45 

-150. 

2 

11. 

38 

-81.  2 

4.  71 

-1.5 

4.64 

0.  2 

6.27 

173. 

1 

10. 

18 

-154.  3 

4.80 

-53.9 

4.  90 

-52.  2 

6.  05 

137. 

6 

8. 

47 

143.  1 

5.  35 

-106.  1 

5.  36 

-103.3 

5.86 

103. 

2 

7. 

13 

89.  7 

5.84 

-152.2 

5.87 

-150.5 

5.51 

70. 

3 

6. 

35 

41.0 

6.25 

159.0 

6.  30 

160.  2 

4.99 

37. 

0 

6. 

21 

-5.  3 

6.99 

112.  1 

6.99 

113.  1 

4.68 

3. 

0 

6. 

51 

-51.8 

7.  49 

66.6 

7.54 

68.4 

4.48 

-30. 

7 

7. 

33 

-102.0 

7.25 

19.4 

7.  32 

22.  0 

4.43 

-65. 

6 

8. 

92 

-163.  3 

6.65 

-26.9 

7.04 

-25.9 

4.94 

-101. 

3 

11. 

29 

112.5 

5.81 

-79.  3 

7.  10 

-73.4 

5.  75 

-133. 

5 

12. 

01 

12.6 

6.31 

-133.  3 

6.46 

-121.  3 

6.  18 

-165. 

5 

11. 

59 

-61.5 

7.25 

179.  3 

7.  30 

174.6 

6.91 

160. 

7 

9. 

15 

-93.  7 

6.90 

135.8 

6.45 

133.  3 

7.83 

125. 

8 

6. 

27 

-142.5 

6.56 

85.4 

6.  35 

83.  2 

9.  07 

87. 

1 

5. 

53 

169.5 

6.  25 

34.5 

6.  00 

34.  1 

1.  78 

41. 

0 

6. 

07 

126.  3 

5.85 

-17.5 

5.61 

-14.6 

12.72 

-22. 

9 

7. 

54 

87.  7 

5.50 

-69.5 

5.  30 

-64.8 

12.53 

-104. 

0 

9. 

91 

51.8 

5.  39 

-121.9 

5.  34 

-118.0 

10.  00 

-174. 

9 

10. 

46 

23.  7 

5.66 

-176. 7 

5.  70 

-171.  1 

7.82 

131. 

2 

11. 

95 

44.  3 

6.05 

131.4 

6.  10 

134.9 

6.53 

87. 

5 
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'jV'iV'iViV'iV'jV'sV'iV'sV'jV'jV'jV'jViV'iViV'jViWfVfjViVVoV'jVVfiViVVfiViV 

THE  INPUT  DATA  FILE  FOR  DISK  10 
3  DIFFERENT  LOAD  FOR  FIXTURE  B 

Vr'iV'jV'jV'iV'j'r-V-'f-ViV's'riVVfiV'iV-V'sV'jV'jVVr'jVVcVrVfVfVfs'fVnVVfsV 


1. 

97 

-117. 

8 

2.  32 

-26. 

9 

1. 

41 

78. 

5 

1. 

89 

-176. 

2 

2.46 

-91. 

9 

1. 

71 

28. 

1 

2. 

00 

130. 

2 

2.44 

-155. 

5 

2. 

16 

-26. 

2 

2. 

27 

81. 

9 

2.51 

147. 

3 

2. 

48 

-84. 

3 

2. 

40 

34. 

8 

2.41 

96. 

6 

2. 

69 

-145. 

9 

9 

62 

-13. 

9 

2.  25 

50. 

1 

2. 

85 

152. 

4 

3. 

04 

-67. 

3 

2.  35 

6. 

9 

2. 

85 

94. 

8 

3. 

87 

-130. 

3 

2.  28 

-35. 

9 

2. 

75 

43. 

2 

5. 

05 

151. 

6 

2.51 

-78. 

5 

2. 

77 

-2. 

9 

5. 

41 

64. 

0 

2.  69 

-123. 

9 

2. 

68 

-44. 

3 

4. 

91 

-11. 

7 

3.  25 

-175. 

1 

2. 

68 

-82. 

4 

3. 

87 

-67. 

6 

4.  27 

126. 

2 

1. 

93 

-121. 

1 

2. 

90 

-118. 

9 

5.24 

56. 

1 

1. 

91 

-162. 

1 

2. 

55 

-166. 

3 

5.  38 

-20. 

6 

2. 

02 

157. 

0 

2. 

68 

147. 

0 

4.49 

-90. 

2 

2. 

42 

113. 

2 

3. 

33 

96. 

3 

3.49 

-147. 

9 

3. 

21 

64. 

6 

4. 

95 

36. 

4 

2.  76 

163. 

0 

4. 

49 

7. 

5 

7. 

33 

-43. 

4 

2.43 

118. 

9 

5. 

70 

-60. 

8 

7. 

51 

-139. 

7 

2.  36 

77. 

5 

5. 

80 

-134. 

9 

5. 

35 

143. 

6 

2.55 

36. 

2 

4. 

89 

158. 

8 

3. 

67 

87. 

9 

3.  32 

-8. 

6 

3. 

72 

105. 

1 

SO 


APPENDIX  F.     FORTRAN  PROGRAM  FOR  DE-EMBEDDING  THE 

FILTER 


DIMENSION  X(21),MDB11(21),PSE11(21) 

MDB12(21),PSE12(21) 

MDB21(21),PSE21(21) 

MDB22(21),PSE22(21) 
REAL  MDB11.MDB 12, MDB21.MDB22 

MAG11,MAG12,MAG21,MAG22 

IM11,IM12,IM21,IM22 
COMPLEX   S11,S12,S21,S22 

A11,A12,A21,A22 

B11,B12,B21,B22,W 

D11,D12,D21,D22 

C11,C12,C21,C22 

F11,F12,F21,F22 

M11,M12,M21,M22 

H11,H12,H21,H22 

PI=4.  *ATAN(  1.  ) 

DO  10  1=1,21 

INPUT  DATA  FILE  IS  DISK8 
(FIXTURE  1  ONLY) 

READ(8,20)DB11,DB12,DB21,DB22,ANG11,ANG12,ANG21,ANG22 
20    FORMAT(8(F7.2)) 

Vll=10.**(-DBll/20.  ) 
V12=10. **( -DB 12/20.  ) 
V21=10.**(-DB21/20.  ) 
V22=10.**(-DB22/20.  ) 
RAD11=PI/180. -ANG11 
RAD12=PI/180.*ANG12 
RAD21=PI/180.*ANG21 
RAD22=PI/180.*ANG22 

QRE11=V11*C0S( RADII) 
QRE12=V12*C0S(RAD12) 
QRE21=V21*COS(RAD21) 
QRE22=V22*COS(RAD22) 
QIM11=V11*SIN(RAD11) 
QIM12=V12*SIN(RAD12) 
QIM21=V21*SIN(RAD21) 
QIM22=V22*SIN(RAD22) 

A11=CMPLX(QRE11,QIM11) 
A12=CMPLX(QRE12,QIM12) 


A21=CMPLX(QRE21,QIM21) 
A22=CMPLX(QR£22,QIM22) 

"S"PARAMETERS  TO  "^'PARAMETERS 

B11=-(A11*A22-A12*A21)/A21 
B12=  A11/A21 
B21=-A22/A21 
B22=l/A21 

INVERSE  OF  FIXTURE  1  "T" PARAMETER 

W=  1/(B11*B22-B21*B12) 

C11=W*B22 

C12=-W*B21 

C21=  -W*B12 

C22=  W*B11 

CALL  OF  'T' PARAMETERS  OF  BOTH  FIXTURE  +  FILTER 

CALL  TMFS(D11,D12,D21,D22,I) 

(INVERSE  OF  FIXTURE  1)  *  (BOTH  FIXTURE+DEVICE) 

H11=C11*D11+C21*D12 
H12=C12*D11+C22*D12 
H21=C11*D21+C21*D22 
H22=C12*D21+C22*D22 

CALL  OF  INVERSE  'T1 PARAMETER  OF  FIXTURE  2 

CALL  FIX2(F11,F12,F21,F22) 

CASCADE  OF  FIXTURE  1,  FILTER  AND  FIXTURE  2 

M11=H11*F11+H21*F12 
M12=H12*F11+H22*F12 
M21=H11*F21+H21*F22 
M22=H12*F21+H22*F22 

"TMPARAMETERS  TO  "Sf 'PARAMETERS 

S11=M12/M22 

S12=  M11-(M12*M21)/M22 

S21=  1/M22 

S22=  -M21/M22 

QIM1=AIMAG(S11) 

QRE1=REAL(S11) 

QIM2=AIMAG(S12) 

QRE2=REAL(S12) 

QIM3=AIMAG(S21) 

QRE3=REAL(S21) 

QIM4=AIMAG(S22) 

QRE4=REAL(S22) 
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C     DE -EMBEDDED  OUTPUT 

MAG11=CABS(S11) 

MAG12=CABS(S12) 

MAG21=CABS(S21) 

MAG22=CABS(S22) 

MDB11(I)=ABS(20*AL0G10(MAG11)) 

MDB12(I)=ABS(20*ALOG10(MAG12)) 

MDB21(I)=ABS(20*ALOG10(MAG21)) 

MD322(I)=ABS(20*ALOG10(MAG22)) 

PSE11(I)=ATAN2(QIM1,QRE1)*180.  /PI 

PSE12(I)=ATAN2(QIM2,QRE2)*180.  /PI 

PSE21(I)=ATAN2(QIM3,QRE3)*180.  /PI 

PSE22(I)=ATAN2(QIM4,QRE4)*180.  /PI 

N=5680 

X(I)=N+320.*FLOAT(I) 

WRITE(6,100)  X(I),MDB11(I),MDB12(I),MDB21(I),MDB22(I), 
1  PSE11(I),PSE12(I),PSE21(I),PSE22(I) 

100   FORMATC  '  ,2X,F6.0,8(1X,F7.2)) 
10    CONTINUE 

END 


Q  iVVriVVfVfVrVrVciViV^ViV^VfiViVVfiViV^ViViV^')V')V,)V'A,'5V'5V,)V,A,')V')'r,iV')V'>V-ft',)V,A' 

SUBROUTINE  TMFS  (Dll ,D12 ,D21 ,D22 , I) 

COMPLEX  A11,A12,A21,A22 
COMPLEX  011,012,021,022 
PI=4.*ATAN(1.  ) 

INPUT  DATA  FILE  IS  DISK  9 
(BOTH  FIXTURE  +  DEVICE  ) 

READ(9,21)DB11,DB12,DB21,DB22,ANG11,ANG12,ANG21,ANG22 
21    FORMATC 8(F7.  2)) 

VI 1=10. **(-DB 11/20. ) 
V12=10.**(-DB12/20.  ) 
V21=10.**(-DB21/20.  ) 
V22=10.**(-DB22/20.  ) 
RAD11=PI/180.*ANG11 
RAD12=PI/180.*ANG12 
RAD21=PI/180.*ANG21 
RAD22=PI/180.*ANG22 

QRE11=V11*C0S( RADII) 
QRE12=V12*C0S(RAD12) 
QRE21=V21*C0S(RAD21) 
QRE22=V22»C0S(RAD22) 
QIM11=V11*SIN(RAD11) 
QIM12=V12*SIN(RAD12) 
QIM21=V21*SIN(RAD21) 
QIM22=V22*SIN(RAD22) 
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A11=CMPLX(QRE11,QIM11) 
A12=CMPLX(QRE12,QIM12) 
A21=CMPLX(QRE21,QIM21) 
A22=CMPLX(QRE22,QIM22) 

D11=-(A11*A22-A12*A21)/A21 

D12=  A11/A21 

D21=-A22/A21 

D22=l/A21 

RETURN 

END 


SUBROUTINE  FIX2(F11 ,F12 ,F21,F22) 

COMPLEX  A11,A12,A21,A22 
COMPLEX  B11,B12,B21,B22,W 
COMPLEX  C11,C12,C21,C22 
COMPLEX  F11,F12,F21,F22 
PI=4.  *ATAN(  1.  ) 

INPUT  DATA  FILE  IS  DISK  10 
(FIXTURE  2  ONLY) 

READ(10,20)DB11)DB12JDB21>DB22)ANG11,ANG12JANG21,ANG22 
20    F0RMAT(8(F7.2)) 

Vll=10.**(-DBll/20.  ) 
V12=10.**(-DB12/20.  ) 
V21=10.**(-DB21/20.  ) 
V22=10.**(-DB22/20. ) 
RAD11=PI/180.*ANG11 
RAD12=PI/180.*ANG12 
RAD21=PI/180.*ANG21 
RAD22=PI/180.*ANG22 

QRE11=V11*C0S( RADII) 
QRE12=V12*COS(RAD12) 
QRE21=V21*COS(RAD21) 
QRE22=V22*COS(RAD22) 
QIM11=V11*SIN(RAD11) 
QIM12=V12*SIN(RAD12) 
QIM21=V21*SIN(RAD21) 
QIM22=V22*SIN(RAD22) 

A11=CMPLX(QRE11,QIM11) 
A12=CMPLX(QRE12,QIM12) 
A2 1=CMPLX( QRE2 1 , QIM2 1 ) 
A22=CMPLX(QRE22,QIM22) 

B11=-(A11*A22-A12*A21)/A21 

B12=  A11/A21 
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B21=-A22/A21 
B22=l/A21 

W=  1/(B11*B22-B21*B12) 
F11=W*B22 

F12=-W*B21 
F21=  -W*B12 

F22=  W*B11 

RETURN 

END 
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